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CHAPTER I 


INTRODUCTION 


The frequency of discovery of gas condensate reservoirs, 
as well as their importance to the industry, has increased 
greatly in recent years with the trend toward drilling into 
the deeper environment of higher temperatures and pressures. 
Craft and Se point out well test data confirm this 
trend as shown in Figure |. This graph shows the discovery 
trend for 17 parishes in southwest Louisiana. The reservoirs 


were separated into oil and gas or gas-condensate types on 


DEPTH, THOUSANDS OF FEET 





NUMBER OF DISCOVERIES 
PER 1000 FT INTERVAL 


Пе си А. Discovery frequency of 
oil and gas or gas-condensate 
reservoirs versus depth. (erat 


and Hawkins, Applied Petroleum 
Engineering.) 
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the basis of well test gas oil ratios and the API gravity of 
the produced liquid. While oil discoveries predominate at 
depths less than 8000 feet, gas and gas-condensate discover- 
ies predominate below 10,000 feet. The decline in discover- 
ies below 12,000 feet is due to the fewer number of wells 
drilled below that depth rather than a drop in the occurrence 
of hydrocarbons, 

Many of these high pressure condensate reservoirs are 
discovered in conditions above their dew points. With pres- 
sure reduction, when placed on production, they allow retro- 
grade condensation of the valuable heavier hydrocarbons within 
the reservoir rock. Once condensation has taken place within 
the reservoir, this liquid cannot be produced without con- 
siderable expense. The mechanisms by which this condensation 
may be retrieved is either by revaporization or by repressuri- 
zation by injected dry gas. However, several pore volumes of 
dry gas must usually be injected before sufficient revaporiza- 
tion occurs, making this process costly. Even then, the dry 
gas may not sweep all the volume in which condensation has 
taken place. In order to retrieve this condensate, which is 
many times more valuable than the gas itself, or to prevent 
its condensation within the reservoir, it is necessary to in- 
vestigate the possibility of dry gas cycling and pressure 
maintenance. 

Cycling of the dry residue gas may act to maintain the 


pressure in the reservoir at or near the dew point, thereby 
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preventing the occurrence of the retrograde condensation 
phenomena. It may also tend to revaporize the condensate and 
sweep it, along with the cycled gas, out of the reservoir. 

Even though it is apparent that cycling a reservoir 
which is susceptible to the retrograde phenomena will produce 
more of the valuable hydrocarbons in place, the economic fac- 
tors involved are all important in deciding whether to invest 
in this procedure. Roughly, a reservoir that contains as much 
as fifty barrels of condensate per million standard cubic feet 
of gas is a good prospect for gas cycling. On the one hand, 
the dry gas sales will be postponed until the cycling is prac- 
tically complete, resulting in decreased present value of the 
gas. However, it is generally possible to produce the con- 
densate faster, hence earlier, when cycling is adopted and 
more dollars are generally obtained due to the increased ulti- 
mate recovery. The purpose of this thesis is to evaluate whether 
or not these advantages are worth the greater investment in- 
volved for the field studied. 

Data are obtained from an actual field on the Gulf Coast 
of Louisiana containing approximately 600 billion SCF of gas. 
The size of the reservoir was determined from geologic data 
and from production data to date. Future prediction of reser- 
voir performance was made for pressure depletion, normal cycling 
and pressure maintenance methods of exploitation. An economic 


comparison of the cases studied is presented. 





CHAPTER II 


DISCUSSION 


A. Nature of Gas Condensate Reservoirs 


Gas reservoirs are defined as having a lean gas, con- 
taining a minimum amount of the heavier hydrocarbons, where 
the composition of this produced gas does not change with 
pressure reduction as the reservoir is depleted. Oil reser- 
voirs range in gas content from zero, or dead oil, to a gas- 
oil ratio of a few thousand cubic feet per barrel. In gas 
reservoirs, on the other hand, the liquid may be vaporized 
into 100,000 cu. ft. per barrel or more. Gas condensate 
reservoirs are classified between these two extremes. Gas 
condensate reservoirs, or distillate reservoirs as they are 
Sometimes calľed, usually produce light colored or colorl ass 
stock tank liquids with gravities above 45 degrees A.P.I. at 
gas-oil ratios in the range of 3,000-100,000 SCF/bbl. The 
production from these reservoirs is predominately gas from 
which some liquid is condensed in the surface separators, 
which phenomenon gives them the name gas-condensate. 

Another definition of the three types of fields described 
can be derived from an interpretation of the pressure-temperature 
phase diagrams. (See Figure 2 for a typical diagram.)  When- 


ever a reservoir fluid is discovered at a temperature above the 





5 
cricondentherm, or maximum two phase temperature for the fluid, 
the composition of the produced fluid will not change as the 
reservoir is depleted since the reservoir temperature can be 


expected to remain fairly constant (Path A-A, on the diagram). 
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Fig. 2. Pressure-temperature phase 
diagram ofa reservoir fluid.: ИСтагв 
and Hawkins, Applied Petroleum Reservoir 
Engineering. ) 


There may, however, be liquid produced as the fluid flows up 
the well bore and into the separators if temperature reduction 
below the circondentherm occurs within the tubing string (Path 


A-A, Om траве 2). If. on the- other hand, the reservon Tiis 





found to be at a temperature below the critical temperature 

(T, ) and pressure above the bubble point, the reservoir would 
be in a one phase liquid state. This is called a bubble 

point reservoir and the fluid will remain a liquid although 

as the bubble point pressure is reached, with pressure deple- 
tion in the reservoir, gas will come out of solution in the 

oil and will be produced along with the oil in an ever increas- 


ing quantities (Path C-C, in Figure 2). In dew point or retro- 


| 
grade gas condensate reservoirs the reservoir fluid exists at 

a temperature between the critical temperature and the cricon- 
dentherm and at a pressure above the dew point pressure of the 
fluid. In this situation, when the pressure declines, composi- 
tion of the reservoir fluid, which may be termed gas in this 
state, will remain the same until the dew point is reached. 

At this point as the fluid expands isothermally, a liquid be- 
gins to condense out of the gas and will adhere to the walls 

of the rock, tubing, etc., leaving a gas leaner in the heavier 
hydrocarbons to be produced, The term retrograde is given to 
this phenomenon since one would expect vaporization to take 
place with a pressure decline rather than condensation. At 
some intermediate pressure, between the dew point and abandon- 
ment, this retrograde phenomenon will reach a maximum, below 
which pressure the condensate will revaporize into the gas. 
(Pressure-Temperature phase relationships for gas-condensate 


. ) 


reservoirs are shown on Figure 2 by the path from B to B, 





B. Retrograde Loss of Condensate 


Liquid that is condensed inside the reservoir rock 
remains immobile in the usual case and therefore is lost to 
production unless or until revaporization can take place 
either by further pressure decline, repressurization, or by 
sweeping with a dry gas. The alternative to revaporization 
of this condensed liquid is to prevent it from condensing in 
the first place and this possibility leads to the investiga- 
tion of pressure maintnenace by injection of dry gas into the 


reservoir. 


с Experience in Gas Cycling 


In general, the gas with the greater condensate content 
exhibits the greater condensate loss upon pressure depletion 
so usually it proves to be the more profitable reservoir to 
cycle with gas. In gas cycling the condensed liquid is sepa- 
rated from the wet (produced) gas in a gasoline plant or other 
surface facilities. The residue or dry gas is then re-injected 
into the reservoir where it maintains the reservoir pressure 
and also serves to drive the wet gas toward the producing wells. 
This procedure is continued until the dilution by dry gas at 
the producing well renders the operation unprofitable at which 
Eime the reservoir is blown down, or allowed to produce nor- 
mally through pressure depletion. During this pressure blow- 


down not only will the lean gas in the reservoir be produced 





but also additional liquid may be recovered from both the 
dry gas invaded zone and the uninvaded portion of the 
reservoir. 

Bee has shown that cycling in the Grapeland Field, 
Houston County, Texas, resulted in a recovery of seventy-five 
per cent of the initial butanes and heavier hydrocarbons in 
the reservoir, while straight pressure depletion would have 
recovered only thirty-nine per cent of the initial butanes-- 
plus. The liquid that could have been recovered by pressure 
depletion, then, was only fifty-two per cent of thatewhieb 
was actually produced by the cycling process, In this example 
it was shown that the ultimate liquid recovery was reached in 
approximately eleven years but twenty-one years would have been 
required by straight depletion, thereby making the present worth 
of the liquid products much greater with gas cycling. 

In a study of the Krotz Springs Field, St. Landry Parish, 
Louisiana, Be تی‎ predicted that twenty per cent more con- 
densate could be recovered by cycling, resulting in a gross 
cash gain after taxes of $2,000,000 for the twenty-five year 
Period studied. 

сай. ےت‎ concluded that benefits in improved 
recovery of stock tank condensate by pressure maintenance 
methods may be as great as 300 per cent in rich condensate 
reservoirs but grade progressively downward to negligible 


benefits for the lean condensate reservoirs with GOR's 


Sreaters than 300,000 cubic foot/bbl. 





Although cycling may appear to be the ideal solution 
to the retrograde condensation problem there are other fac- 
tors that may make this procedure uneconomical. The deferred 
income from the dry gas sales may prove to be a significant 
factor in the economic analysis of the field. Another eco- 
nomic consideration is the added investment required for 
additional injection wells, gas compression and distribution 
system to return residue gas to the wells and for a liquid 
recovery plant. Also, it must be remembered that even with 
gas injection at a pressure above the dew point all of the 
liquid hydrocarbons may not be recovered. 

Three different recovery factors must be applied with 
gas cycling. The microscopic displacement efficiency is 
approximately 70 to 90%. The volumetric sweep efficiency, 
or the per cent of initial pore volume invaded by the sweep- 
ing gas to abandonment of the producing wells can vary from 
БОВ Со 905. Finally a permeability stratification factor 
must be applied to take care of the problem of gas sweeping 
through the more permeable strata and reaching the producing 
wells before the tighter stringers can be swept. The over- 
all condensate recovery factor for a gas cycling project is 
the product of these three factors. On the plus side of 
recovery by cycling is the fact that additional liquid is 
produced during blow-down from the less permeable or unswept 


portions of the reservoir as these products bleed into the 
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more permeable streaks and ultimately to the wells. Also the 
efficiency of extracting condensate from the wet gas is higher 
with the type of separator plant that is generally used ina 
cycling project. 

Any reservoir exhibiting retrograde properties and hav- 
ing a richness of 50 bbls/MMSCF of condensate is a good pros- 
pect for cycling of the dry residue gas to improve recovery. 
Each reservoir, however, is an individual case and, as such, 
requires a separate study. Variables such as permeability 
stratification, reservoir geometry, gas properties, and market 
availability can team up to make one reservoir an extremely 
profitable enterprise while another may prove entirely infea- 
sible. This thesis presents the study of one field based on 


the fluid properties and geologic structure of that one field. 





CHAPTER III 


DEFINITION OF THE FIELD STUDIED 


A. Development of the Field 


Production from this field began in January, 1961 and 
from that time two to three wells have produced continuously 
as shown in Table II. The pool is one in which higher than 
normal pressures are encountered. It has been obvious from 
the low gas-oil ratios measured that the produced gas was 
rich in condensate. By January 1, 1964, a total of 20 billion 
standard cubic feet of wet gas had been produced from the two 
wells completed in the zone studied. The pressure had dropped 
from an estimated initial pressure of 8838 psia to 8100 psia 
as shown in Figure 7. Based on these production data, calcu- 
lation of the initial gas in place yielded the value of 606.9 
billion standard cubic feet. These values form the basis for 


future performance predictions from the base date of January, 


1964. 
B, Structure and Reserves 

Interpretation of available geologic data indicates a 
deep-seated salt domal structure. Producing depths at which 


this study was conducted were from approximately 11,300 to 
11,600 feet. A structure contour map of the top of the pro- 


ducing formation is presented as Figure 3. The productive 


1 
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13 
interval is taken to be a shaly-sand with intergranular 
porosity. 

There are a total of nine wells drilled through the 
formation studied and from which data are available from 
various kinds of logs. Porosity was determined using micro- 
logs from some wells and sonic logs from others. The porosity 
throughout the formation is fairly constant, varying from 257 
to 29% in the pay zone. Although shaly streaks were noted, 
this fact was taken into consideration in counting sand, 
wherein the shale sections were discounted. For purposes of 
the material balance calculations an average uniform porosity 
of 28% was used. 

Electric logs or induction-electric logs were available 
from all nine wells. Three of these logs are shown in Figure 4. 
These logs were used to compute formation water resistivity 
using the Spontaneous Potential Curve. Water saturation was 
calculated from the induction log or lateral curve as avail- 
able. Water resistivity was found to be approximately .03 ohm- 
meters at formation temperature and this value was used through- 
out the field. Computed water saturation tended to vary from 
8-10% near the top of the gas zone to about 25% nearer the 
bottom of the pay and it grades to 100% water saturation in 
those wells where a gas-water contact is observed. A value of 
20% water-saturation was used for the entire formation as being 


the weighted average value. 
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Side wall cores were available from one well. The 
porosity values obtained from the sidewall cores checked pre- 
cisely the values obtained from the sonic logs and micrologs. 
The permeability, however, did not appear uniform, nor did 
these values check well against the effective permeability 
calculated from a pressure draw-down test or against permea- 
bility calculated from a resistivity gradient observed in some 
of the wells. The six values obtained from side-wall cores 
were 84, 134, 50, 512, 131 and 116 millidarcies. The effec- 
tive permeability to gas obtained in the pressure draw-down 
test proved to be 50 millidarcies and this was the value used 
in the material balance calculations since the method of tak- 
ing side-wall cores tends to result in fractured, non- 
representative samples of the reservoir. A plot of the data 
taken from the pressure draw-down test is shown in Figure 5. 

For volumetric sweep-out calculations, 50 millidarcies 
was used for permeability while the effective thickness of pay 
was allowed to vary over the field. This gave a value of flow 
capacity for each point in the reservoir which was used in the 
sweep efficiency program. 

Connate water saturation was found to be 20 per cent. 
This was an average value based on water saturation in the 
lower portion of the gas zone. Residual gas saturation to 
water is estimated to be 30 per cent. This is justified by 


the gas saturation near the top of the water column which, it 
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12 
can be assumed, has encroached into the reservoir to some 
extent during a recent producing time. 

Determination of net effective thickness in this forma- 
tion was very difficult because of the shalyness of the sand 
and because of the meager geologic information. Only nine 
wells have penetrated the 2000 acres of estimated pay forma- 
tion. The locations of the two faults--one thought to have a 
closure of about 100 feet in the northern extremity of the 
pool, the other with an estimated closure of about 120 feet 
in the south--are by no means exactly defined. An error of 
the position of the fault to the north of only a few feet 
would produce a substantial error in the estimated reserves 
because of the great thickness of pay adjacent to this fault. 
Calculation of reserves based on the volumetric method is 
presented in.Table II. The east and west boundaries of the 
field are located at a gas-water contact which has shown 
little or no encroachment during the current gas production. 

An isopach map produced from the nine control wells is 
presented as Figure 6. 

A volumetric gas balance was made using the production 
data to date and bottom hole pressure tests taken in 1961 and 
1963. A bottom hole pressure was obtained from the draw-down 
and build-up pressure test conducted in connection with this 
study. This gas balance was calculated on a Control Data 
Corporation 1604 high speed computer using a program developed 


by the Union Oil Company of California. 
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19 
J. Donald Clark, Chief Reservoir Engineer for the Gulf 
Division of the Union Oil Company of California, has developed 
and applied or supervised the application of a very simple 
transient energy balance equation to numerous natural water 
drive gas reservoirs found in the Gulf Coast arae A primary 
requirement of this equation is the derivation of the gas for- 
mation volume factor, БЕ to convert standard cubic feet of 
gas to reservoir barrels. This conversion for any particular 


reservoir is presented mathematically as follows: 


MEM 800 08 +٦ 


(٦‏ ۹+ 9و9 


Where: 15.025 is base pressure for Louisiana (psi) 
520 is standard temperature OR 

.61 converts cubic feet to barrels 

is reservoir temperature, degrees Rankine 

is gas deviation or compressibility factor 

is the reservoir pressure, psia 

as a subscript will denote initial or 
original conditions 


к> ку HU 


The cumulative reservoir barrels of gas produced or withdrawn 
is equal to the cumulative reservoir barrels of gas expansion 


in the reservoir, or 
2. GG В = G(B -B ,) 
pw 5 5 81 
Where: Gpy is cumulative standard cubic feet of 
wet gas produced 


G is standard cubic feet of original gas 
IG place 


In the case of water drive reservoirs and water production 


the statement of the equation then becomes: (Cumulative Net 
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Reservoir Withdrawals) - (Cumulative Reservoir Expansion) = 


(Cumulative Water Influx), or 


3. (G B -*W)- G(B - B,) -W 
с. ©. ) 


gi e 
4. G B + W - G(B - B_.) + АХАР 
( ae 3) ( 5 T 0. 
The water influx equation, ШЕ = ALAPQ,, is a simplified 


method of writing the equation presented by gen in his 
paper entitled "Water Influx into a Reservoir and Its Appli- 
cation to the Equation of Volumetric Balance." Details of 
solving the water influx equation were presented by van 
Everdingen and SON The "A" is the water influx constant 
which encompasses conversion factors, porosity, radius of 
reservoir, effective reservoir thickness, effective compressi- 
٣7٢٦ہ‎ reservoir fluid, viscosity of reservoir fluid ang 
fraction of reservoir perimeter exposed to water drive. The 
above parameters are constant for any particular reservoir 
and when used with the volumetric balance equation can be 
solved as an all encompassing constant for that particular 
reservoir. 

The true simplification is in the manner of arriving at 
Ehe value of ZAPQ, for various balancing periods. 

The dimensionless time equation, En = KT/$u0 x, is 
derived from the hydraulic diffusivity factor, а - K/yugC 


presented by An in 1933. The K, or permeability factor, 


Becually is the one moist difficult tos/measure since it is 
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effective permeability to water, which we seldom if ever 
obtain. This value could be found by running draw-down tests 
in salt water wells, but actually it can be estimated to a 
Sufficient degree of accuracy. Viscosity of water can be 
obtained from various sources to a reasonable order of magni- 
erde, Therefore, it becomes important only to establish an 


order of magnitude value of t in order to extract reasonable 


D’ 
values of the corresponding в. the dimensionless water influx, 
from the tables. 

Experience has shown that a proper time period should 
not exceed three months or approximately a 90-day producing 
period when attempting to strike a balance, especially if one 
plans to calculate future predictions. The constant "A" be- 
comes an over-all correcting constant of transient flow, and 
allows us to solve very simply a satisfactory value of original 


hydrocarbons in-place and the total water influx at any time in 


the producing life of the reservoir. 


Cumulative Production and Pressure History 

Production from the reservoir is tabulated on a monthly 
rate and cumulated on a three-month basis. These data are 
broken down between oil, gas, condensate and water (both 
fresh and salt water). When rates of production have consider- 
meee variation, a monthly production rate should be plotted as 
shown on Figure /7. This plot is important im order to show 


the effect of rate on pressure when natural water drives exist. 
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Figure 7 
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Bottom hole pressure measurements are plotted versus 
time on the same graph as production rate. The length of 
shut-in time prior to each pressure measurement should be 
approximately equal and the individual pressures should have 
reasonably close values. If not, build-up or drawdown tests 
should be run to determine the cause of pressure variation. 
A pressure curve is then constructed through these evaluated 
pressure points. Experience and ingenuity are of tremendous 
value in constructing a reasonable curve. The production 
rate graph must be studied in order to prepare the most logi- 
cal pressure-time curve. Pressure for uniform periods of pro- 


duction are taken directly from this graph. 


Water Influx Data and Basic Calculations 
The next step requires the gathering of basic informa- 
tion to approximate dimensionless time data and corresponding 


dimensionless rates. 


KT(86,400) 
ETE TUR 2 
u$(0.5 x 107^) р, (930) 


Here geological data can be used to help determine the radius 


of the hydrocarbon reservoir. Only a close approximation is 
necessary. Compressibility of salt water can be found in the 
literature. Viscosity of salt water can be approximated by 


using chemical handbooks, available texts, and various publi- 


cations; porosity is approximated from core analyses or logs 
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and permeability to water can be estimated from air permea- 
bility of cores analyzed. 

Dimensionless time values are calculated for cumulative 
three-month periods. The corresponding са values are taken 
from the table. Integration of the values for AP times Q. 
over the cumulative time intervals yields the value which, 
when multiplied by the constant A, provides the cumulative 


water influx. 


Gas Conversion to Reservoir Barrels 

The vapor equivalent of the distillate and fresh water 
produced must be added to the amount of dry gas produced. The 
constant for converting liquid water to the vapor state is 
7368.48 cubic feet per barrel. The distillate conversion fac- 
tor was calculated from the hydrocarbon analysis of the sepa- 
rator liquid shown on page 89, Appendix II. This factor was 
found to be 768 standard cubic feet of vapor per barrel of 
distillate. 
Method of Least Mean Squares for Striking an 
Energy Balance 

The method of least mean squares is used for striking 


an energy balance for the natural water drive gas reservoir, 
m + Ше - е - D = A (ZAPQ, ) 


The two unknowns are G, the original wet gas in-place, 


апа A, the water influx constant. The method herein presented 





25 
requires the solution for single values of A and G to meet all 
periods of past producing time. We, therefore, can arrive at 
two equations with two unknowns and place proportional weight 


to the length of production time: 


2 
[7 (6, + و‎ - GIG, - Bi) (ZAPQ,)] = AA | 


w 


2 
E. T М 6) (В - B .)] - C[ @ - B .) 


g gi g si)] = A[(ZzAPQ )(B -B ;)] 


gi 

These two equations can then be solved for A, the water 
influx constant, and G, the original wet gas in-place. 

In order to check the validity of the values of A and 
G obtained by this method, multiply the LMSA value of A times 
LAPQ, and solve for the apparent value of G. This is called 
"periodic check of balance for G assuming LMSA values of A 
as constant." The van Everdingen and Hurst solution to the 
diffusivity equation is here applied to single-phase fluid 
transient flow into a well bore for determination of the rate 
and volume of encroachment of water. The gas field is assumed 
to be embedded into an aquifer of large extent. In this 
instance, the reservoir is taken to be a large well, the 
radius of which is the equivalent field radius, and the aqui- 
fer is considered to be the reservoir from which fluid is 
withdrawn and water is oe 

The computer program, indicating method of computation, 


Weservoir conditions, initial gas in place, and cumulative 
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water influx from past production data, is included as 
Appendix I. 

Since the available geologic information is considered 
isufficient to describe the reservoir accurately enough for 
calculation of reserves, the results obtained from the volu- 
metric gas balance were used for estimating original gas in 
place and water encroachment. These results were also the 
basis for the future performance prediction and economic 
analysis of alternate methods of exploitation. The original 
gas in place was found to be 606.9 billion cubic feet measured 
at standard conditions (15.025 psi at 60°F). Water influx was 


found to be negligible. 


Cc. Reservoir Fluid Characteristics 

Samples of both the separator liquid and gas were ob- 
tained in March, 1964, and recombined in accordance with the 
measured gas-liquid ratio (stock-tank) of 17,150 SCF separator 
gas per NM ck tank liquid. It was difficult to obtain 
constant gas-oil ratios for the purpose of getting a repre- 
sentative sample and flowing at a constant rate for about three 
days was required before a steady GOR was achieved. Analysis 
of the recombined sample indicates a dew point pressure above 
the present bottom hole pressure. Retrograde condensation 
was taking place and two phase flow occurring near the well 


bore caused the fluctuating producing gas-liquid ratio. Reser- 


voir temperature at time of sampling was 240 degrees farenheit. 
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Primary separation temperature was 125 degrees farenheit and 
primary separator pressure was 1020 psig. Dew point pressure 
was found to be 8260 psia and gas richness was 66.07 bbls. of 
stock tank liquid per million standard cubic feet of well 
stream effluent. Detailed results of the P-V-T investigation 
are included as Appendix II. 

Samples of separator gas and separator liquid were ana- 
lyzed by the conventional low temperature distillation method. 
The results are reported shoving both the composition of each 
sample and the computed analysis of the well stream based on 
the gas-liquid ratio of the primary separator. The separator 
liquid production was calculated from the measured stock tank 
production by applying the determined shrinkage factor. 

Following the compositional analyses, portions of the 
primary separator liquid and gas were physically recombined 
in their produced ratio in a variable volume, glass-windowed 
equilibrium cell. Determinations on this mixture were divided 
into the following two main categories: 

J. Dew point pressure determination and pressure- 
volume relations on a constant weight of reservoir fluid 
at the reservoir temperature: The procedure consisted of 
establishing equilibrium between gas and liquid phases at 
a low pressure and measuring the volumes of liquid and 
gas in equilibrium at that pressure. The pressure was 


then raised by the injection of mercury into the cell and 
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phase equilibrium established again at the higher pressure. 
This procedure was repeated until all of the liquid phase 
had vaporized, at which point the saturation pressure was 
observed. The cell pressure was then raised above the dew 
point pressure in order to determine the supercompressi- 
bility characteristics of the single phase vapor. As a 
check on all readings and particularly to verify the dew 
point, the cell pressure was incrementally reduced, equi- 
librium established and volumetric readings made. Reported 
in Appendix II are the relative volume relations and spe- 
cific volumes of the reservoir fluid over a wide range of 
pressures as well as deviation factors of the single phase 
vapor above the dew point, (Also reported in Appendix II 
are the dew point pressures resulting from recombinations 
at gas-liquid ratios above and below the ratio measured 
at the time of sampling.) 

II. Compositions of the produced well stream and the 
amount of retrograde condensation resulting from a stepwise 
differential depletion: This procedure consisted of a 
series of constant composition expansions and constant vol- 
ume displacements with each displacement being terminated 
at the original cell volume. The gas removed during the 
constant pressure displacement was charged to the low 


temperature fractional distillation equipment for analysis. 
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D. Conformance 
The description of the simultaneous flow of fluids 

through porous media in terms of relative permeability and 
capillary pressure has been adequately discussed in the lit- 
erature. For example, Muskat” in Chapter VII briefly dis- 
cussed the dynamical effects associated with capillary 
phenomena. He also outlined the theory of potentiometric 
models and illustrated the application of the potentiometric 
model to cycling systems. Douglas, Peaceman, and and, 
presented a method for calculating multi-dimensional immiscible 
displacement based on a ime i solution of a finite differ: 
ence analogue describing the displacement process. This tech- 
nique insures that at any point in the sand, the flow of oil 
per unit area across the direction of flow can be represented 
by a vector ы. which is given by 


=” ro 
O 
Ho 





(Vp. + 7 g Vh). 
Ih t equals Т ° + 7 ° a Lt 
e operator, y, equ iot j OV ; ‚ ] are uni 
vectors in the x and y directions, respectively; k is the 
local permeability, К о is the relative permeability to oil, 


u the oil viscosity, Ро the pressure in the oil phase, Ро the 





oil density, and g the acceleration of gravity. Similarly, 
> ча, 
= - + 
Ч i. (ур, + P, 8 Vh) 


w 


Fluid densities are assumed constant and oil and water 
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potentials are defined 
۶ ۶ ى‎ 97٤7 


0 7p ІН Ы СЕП 


The author proceeds to solve the problem by a system of dif- 
ference equations. 

In solving this reservoir a technique demonstrated by 
Dr. S. J. Pirson at The University of Texas is followed. In 
this method the reservoir is divided into a grid and each 
intersection of the grid is treated as a point (I, J) within 
the reservoir where Darcy's law may be applied. The flow 
across the grid between these points, q, is now analyzed. 
The rate of flow over a distince Ax across an area open to 


flow, Ay, can be determined by Darcy's law: 





Р - P 
= سج سا‎ 2 
41 70 7 5 چرم‎ а 


In order to satisfy the continuity equation the vector 


sum of the flow in all four directions must equal zero unless 
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there is fluid withdrawal or fluid injection at the point 


(Р. . -P, .)+K 
O av av — ۹93 کت‎ 


(P ) 


VOD и 


+K, ( P 


av 


қР” =P. 
ll To 


It ís seen that without injection or production from the 
point studied Q=0 and the Khay is the average flow potential 
between adjacent points. This equation is now solved for P, 


2 


and we find: 


206۰ 

Е + 

ПИРА МА ЕК л Ы 
و‎ . МИ С m uu 

р Е ЕРТЕ ۷۷ت‎ СД јеле 
i,j 
E 4 n n 
rca c RN LM NE S or 


The boundary conditions are those restrictions on the 
limits of the connected region and the position and magnitude 
of the imposed flows or pressures. At the boundary of the 
permeable sand, the velocity of flow perpendicular to the 


boundary is zero, so it is sufficient to establish equal poten- 


tials across the boundary. Since (VÖ) = O across the boundary, 
NS Flow occurs. This same condition occurs at a fault bound- 
ary. As the water saturation increases near the water table, 


the relative permeability to water increases greatly while 


that to oil diminishes. By assigning these relatively high 
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values of flow potential in the equation aquifer pressure is 
maintained in the advancing flood front and movement of the 
front is regulated across this boundary. 

In solving this potential equation for an oil reservoir, 
the term ا‎ is treated as a constant. However, 
since the flow equation depends on pressure, the constant for 


a gas field is augmented by a variable pressure term: 


(а„)(Р_)(Т)(=) 


K dP 
0111ی‎ 1 - 2٣٦ 
For Oil For Gas 


Volumetric sweep efficiency was investigated by this 
technique with the aid of a mathematical model which was run 
on the CDC 1604 high speed computer. This’ computer program 
is included as Appendix III. The program assumes steady-state 
flow at a point, and sequentially at each point in the reser- 
voir in order that Darcy's flow equation may be applied. Using 
an alternating direction iteration procedure, the computer then 
reduces the total flow at each point in four directions to zero 
and by many iterations establishes each potential from the 
adjacent four potentials in order to achieve a balance. At 
each point selected on the model the quantity of flow into or 
out of the reservoir is assigned in the event a well is located 
aE that point. The flow potential, or the product of effective 
permeability times the thickness of effective pay, is also 


assigned to each point. By applying Darcy's law in each of 





ee ~  —w ж‏ ہے 
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four directions, by summation of these directional flow rates 
to zero, and by solving this equation for the potential at 
each point, a potential distribution map of the reservoir is 
produced. Isopotential contour lines are now drawn on a scale 
map of the reservoir. The primary criterion for equivalence 
between the model and the actual reservoir is the creation of 
geometrically similar potential fields which are determined 
only by the thickness, permeability, and boundary conditions. 
Flow lines were superimposed on the potential map and the 
method of squares utilized by ED is used to define the 
flow tubes through which a constant quantity of reservoir 
fluid would be expected to flow. 

Flow lines of a fluid in a porous medium, as defined 
in Darcy's law, are everywhere orthogonal to the isopotential 
lines. The method of squares employs the fact that if a family 
of curves is constructed such that they are orthogonal to one 
member of another family of curves, then the second family of 
curves may be drawn orthogonal to the first. This is done by 
constructing them such that adjacent distances between the 
streamlines are equal to adjacent distances between the iso- 
potentials. A stream tube is a two-dimensional flow system 
defined as the area between two consecutive stream lines. Flow 


rate can be defined by the following: 


_ k Ad 
Q = (h)(- 1) (Ax) RE 
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In any one square Ax = As and Q is constant so: 


At =# (constant) پر‎ Ger 
k h A 


The increment of time can now be calculated for each of the 
Squares within each flow tube so that the position of the 
front in each flow tube can be drawn at any selected time 
period. This procedure is repeated for various injection- 
production well combinations until the maximum volume of the 
reservoir is swept. Rearrangement of the well positions in 
the reservoir matrix is achieved simply by changing the quan- 
tity of flow into or out of that point in the reservoir. Sweep- 
out patterns were predicted with one, two and three injection 
wells and by using four producing wells. The locations of 
existing wells were used as much as possible but alternate 
new locations were also tried. The optimum combination was 
found to be two injection wells and four producing wells. 

Wells marked "Well 4" and "Well 10" on Figure 8 are to 
be recompleted as injection wells. Two new wells are required 
dtithe locations marked "Well 5" and "Well 11." An existing 
well at location "Well 9" is to be recompleted in the pay zone 
studied since it is now producing from a deeper zone. "Well 6" 
is producing from this formation and will continue to be used 
mene this project. 

A machine print-out of the potentials has been reduced 
to a map shown in Figure 8 picturing the isopotential contour 


lines and on which the frontal position has been drawn at one 
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time just after breakthrough of the injected residue gas. The 
volume of the reservoir swept at this breakthrough time was 
Hound to be 65%. 

The ability of the dry absorber gas to vaporize the 
condensate in the sand was determined by pressure-volume- 
temperature tests run on a recombined sample of the gas. 
Results of these tests are given in Appendix II. The overall 
displacement efficiency of the residue gas was found to be 
02+ 

The effect of bypassing on a local scale, due to permea- 
bility variations, was analyzed as demonstrated by ET 
by arriving at a permeability variation for the sand. Although 
insufficient data are available to establish firmly the perme- 
ability variations throughout the reservoir, a variation .30 
was assigned. This value gives a sweep efficiency at break- 
ط ع۲۳0‎ of 87%. 

Using the product of these three conformance factors, 
the overall volume of the reservoir swept out at breakthrough 
was 51%, with the dry gas cut after breakthrough increasing in 
the manner calculated by ہت‎ in the paper mentioned before. 
Cycling of gas was discontinued when proportion of reservoir gas 
produced was reduced to 67% by dilution with cycled residue gas, 


ПСО the cycling economic limit having been reached. 














CHAPTER IV 


PERFORMANCE PREDICTION 


A. Pressure Depletion Method 


Past production has averaged 25 million standard cubic 
feet per day from the reservoir studied. Until this time no 
firm estimate had been made on the reserves nor had any analy- 
sis been made to predict the future performance and probable 
life of the reservoir. From the past production figures the 
initial gas in place was estimated and the pressure decline 
and water influx was analyzed. From these data it is possible 
to predict the future performance of the field. This was done 
through the use of the computer program included in Appendix IV. 

Because of insufficient pressure decline and the lack of 
back pressure data it is not possible to predict deliverability 
deterioration with declining pressure. It is anticipated that 
as liquid saturation around the well bore increases with static 
pressure reduction and with more severe pressure drawdown, 
deterioration of deliverability will result. For the purposes 
of this study, in absence of definitive data, the rate of pro- 
duction of dry gas was assumed to remain constant and at the 
current rate. The water influx constant obtained from the 
volumetric gas balance program was used to obtain future per- 


formance prediction as well as the calculated quantity of 
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initial gas in place and the adjusted pressure decline 
figures. 

The amount of condensate produced, by year, was calcu- 
lated from the P-V-T data obtained in the material balance 
calculations performed on the recombined gas sample placed 
in the equilibrium test cell. In this cell the pressure was 
depleted by removing successively small quantities of gas from 
the cell. Thermodynamic equilibrium was established, material 
balance calculations were made, and retrograde liquid volumes 
and the gravity of the effluent gas were measured at each 
step. The gas removed in each depletion step was charged to 
the low temperature distillation equipment for analysis. The 
produced volumes of reservoir gas, the determined compositions, 
computed gallons per thousand standard cubic feet (GPM) con- 
tent and respective deviation factors are presented in Appen- 
dix II. The volume of retrograde liquid resulting from this 
depletion study is also shown in Appendix II, both in terms 
of barrels of reservoir liquid and per cent of hydrocarbon 
pore space. Shown also are the compositions of the gas and 
liquid remaining in the reservoir after depletion to abandon- 
ment pressure. The amount of condensate calculated to be pro- 
duced with the well stream at each pressure decrement in the 
future prediction was tabulated. The value of gas sold under 
the terms of the present contract is $.20 per MCF and conden- 


sate is sold for $2.95 per barrel. These prices are assumed 
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Ере гіт for future sales. Operating costs are averaging 
$.00423 per MCF produced and product transportation charges 
for the condensate sold are averaging $.1375 per barrel. 
Application of these values to depletion of the reservoir 
over a fifty-year period with net income discounted at 6% 
results in a present worth of $38,749,249. It is seen that 
in this case of pressure depletion to 1500 psi only 73.51 per 
cent of the initial gas in place was produced. Since the gas 
is to be delivered to a pipeline at a pressure in excess of 
1000 psi, further recovery would require recompression to the 
sales gas line. This secondary recovery process was not con- 
sidered in this study because the value of gas recovered 50-odd 
years hence is not relevant for a study of possible investments 
now. 

It is significant to note that the condensate recovered 
by pressure depletion amounts to 15.4 million barrels which 
is scarcely 38% of the more than 40 million barrels of recov- 


erable oil initially in place. 


В. Pressure Maintenance through Gas Cycling 


A program of gas cycling was considered wherein gas 
sales from the reservoir would be terminated and all gas pro- 
duced from the reservoir, together with enough make-up gas to 
allow for compression and other gas losses, would be reinjected 
after extraction of all condensable heavier hydrocarbons. As 


determined by the sweep efficiency model studies, two wells 
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should be recompleted as injection wells and four wells would 
be used as producing wells as shown in Figure 8. Since the 
installation is located off shore the wet gas would be pro- 
cessed in low temperature units on platforms near the produc- 
ing wells. The residue gas would be piped from the separators 
to the injection platform which would be located at one of the 
injection wells. Double extra heavy pipe would be necessary 
to carry the compressed injection gas to the other injection 
wells since injection pressures in excess of 8000 psi will be 
encountered. 

Condensate produced in the separator units would be 
metered and then reinjected into the gas sales line of the 
gas company where it would be transported with other wet gas 
to the on-shore extraction plant for processing. 

Under pressure maintenance it is assumed that the reser- 
voir will continue to produce a gas of the richness measured 
in the most recent flow tests until breakthrough of the in- 
jected dry gas into the first producing well. From this point 
in time the produced gas would decline in richness in direct 
proportion to the residue gas cut produced. When the amount 
of injected dry gas obtained through the production wells 
reached 1/3 of the total gas processed it was found uneconomical 
to continue gas injection. Injection would be discontinued and 
gas sales would begin to deplete the reservoir to about 


1300 psi or until compression was necessary to charge the sales 
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gas line. Production past this point was not considered part 
о 5 5 Еџау. 

It was assumed that the wet gas produced would be in 
equilibrium at the pressure and temperature produced in the 
Same proportion as that indicated by the equilibrium cell 
tests. In the equilibrium cell the recombined gas sample was 
allowed to expand at constant composition with each pressure 
decrement and then gas was extracted at constant pressure, 
thus simulating reservoir withdrawals. The material balance 
calculations resulting from these tests are believed to depict 
more realistically the actual reservoir performance than calcu- 
lation of the material balance by means of equilibrium constants. 

Based on a permeability variation of .30 as defined by 
В па: те" the dry gas cut at the producing wells was estimated 
from the first breakthrough at a producing well until the eco- 
nomic limit of injection was reached. 

Extraneous gas is available in excess of present avail- 
able market in a reservoir below the zone studied and conse- 
quantly at a higher pressure. This gas would be used as make 
up gas in the pressure maintenance case to compensate for 
losses due to fuel consumption and condensate production. 

After the economic limit is reached it is assumed that 
gas sales will begin again at the present rate and that pre- 
vailing prices would still be applicable. 

Since the field is just beginning to enter the retro- 


grade condensation phase considerably more condensate will be 
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produced by pressure maintenance if this procedure can be 
started at once. The greater the decline in pressure before 
injection commences the greater will be the losses in the 
reservoir and the higher will be the value of extraneous gas 
injected. 

Because of the loss in gas sales during the period of 
pressure maintenance this procedure does not appear economi- 
cally attractive unless market availability for present gas 
sales can be filled from another source of gas. The ultimate 
recovery of condensate through pressure maintenance is excel- 
lent from a conservation point of view since over 90% of the 
available condensable hydrocarbon is expected to be recovered 
through the combination of dry gas injection to maintain pres- 


sure, to be followed by reservoir blowdown to depletion. 


C. Cycling with Pressure Depletion 


Return to the reservoir of available gas from the high 
pressure separator constitutes normal cycling. Since a net 
voidage situation exists as a result of shrinkage and fuel 
losses from the produced gas, normal cycling results in only 
partial pressure maintenance. However, the pressure reduc- 
tion would be little more than 40% during the full economic 
cycling life. 

As in the case studied for full pressure maintenance, 
two wells would be used for injection of residue gas and four 
wells would be used for production of the wet gas effluent as 


indicated in the sweep model study. (See Figure 8.) 
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In this analysis gas sales are continued at the present 
rate which is dictated by the available gas market. Provision 
is made for compressors adequate to compress 40 million stand- 
ard cubic feet of residue gas per day for injection. All wet 
gas produced is processed through the low temperature units on 
offshore platforms. The sales gas and all the condensate pro- 
duced is charged to the sales gas line provided by the gas 
company. 

The amount of condensate produced will be governed by 
the equilibrium saturation at the prevailing producing bottom 
hole pressure until breakthrough of detectable dry gas. From 
this point onward the increasing dry gas cut will be governed 
by the permeability stratification as well as by the areal 
sweep pattern. As the cycling reaches its economic limit, the 
compressors will be shut down and sales gas production will 
continue to deplete the field. Production below 1300 psi was 
not considered in this study, since compression would be re- 
quired to charge the sales gas line which is maintained in 
excess of 1000 psi. 

It is predicted that normal cycling to .67 reservoir- 
gas cut followed by pressure blowdown will result in the recov- 
ery of 26,876,000 stock tank barrels of condensate or 67% of 
the available condensable hydrocarbons. Although the total 
amount of condensate produced through normal cycling is con- 


siderably less than for the pressure maintenance case, the 
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present worth of the income from early gas sales tended to 
make normal cycling more attractive economically. The ad- 
vantage of normal gas cycling, then, is continued income from 
the available market for gas while producing more condensate 
than would have been produced through normal pressure deple- 
tion and by producing this condensate earlier in the life of 


the field and in excess of available market for residue gas. 





CHAPTER V 


CONCLUSIONS 


A. Comparison of Ultimate Recoveries 


In a condensate reservoir, income is derived from the 
sales, not only of the residual gas produced, but also the 
liquids produced from the reservoir gas by mechanical sepa- 
ration or other type of plant processing. The ultimate 
quantity of these fluids produced is of importance as well 
as the rate of this production. 

The original gas in place is the same in any case, 

The ultimate produced quantity of this gas, however, is 
governed by the amount of fuel required in its production, 

or in the production of its condensate, and also by the reser- 
voir pressure at abandonment. The time at which this gas will 
be produced is dependent, of course, on the rate at which it 
is produced which in turn is dependent on the available mar- 
keting possibilities. In each of the alternatives studied the 
ultimate quantity of gas recovered was essentially the same, 
while the timing of the recovery of the gas was dependent on 
the method of field exploitation. 

The amount of condensate recovered is dependent on the 
pressure within the reservoir while being produced, since the 
retrograde properties of the reservoir fluid causes much of 
the heavier liquid hydrocarbon to en within the reservoir 
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rock.: By producing the reservoir at a pressure at which retro- 
grade condensation will not occur or where it approaches some 
minimum value the quantity of condensate ultimately recovered 
is increased. 

The rate at which condensate is produced is dependent 
upon the rate of gas production unless some extraordinary 
means is provided to extract the condensate in advance of gas 
sales, This is the case in most gas reservoirs since pressure 
depletion without gas cycling is generally the method used for 
recovery. 

By cycling a part or all of the gas produced back into 
the reservoir not only can the condensate be extracted in 
advance of a gas market, but also the pressure can be main- 
tained at a level where little or no retrograde loss will 
occur in the reservoir. In the pressure maintenance case 
the rate of condensate recovery is limited only by the size 
of compressor plant installed. In the normal cycling case 
the rate is determined by the sum of the available gas market 
coupled with the capacity of the compressors. Under pressure 
maintenance essentially all, 95%, of the condensate is ulti- 
mately recovered and through normal cycling 67% is recovered. 
When compared with recovery of condensate by pressure deple- 
tion of only 38% of the liquid, it is apparent that cycling 
ís an important conservation measure. 

A comparison of ultimate recovery for both gas and 


condensate is presented in Table I. A plot of predicted 
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cumulative condensate recovery versus time is presented as 


Figure 9. 
B. Economic Analysis of Alternatives 


Since the quantity of gas ultimately to be recovered is 
essentially the same by any method of exploitation, the only 
factor in gas production contributing to its present value, 
and hence important in an economic analysis, is the time at 
which this gas is produced, In order to obtain maximum value 
from this valuable resource, the gas should be produced in as 
great a quantity as can be marketed. In the case of pressure 
depletion and in the case for normal gas cycling, gas was sold 
at the maximum rate Limited only by the market obtainable. In 
the case for pressure maintenance the sale of gas was postponed 
for thirty years and this fact contributed to the reduced pres- 
ent worth for this alternative. 

Through pressure depletion 15,430,000 barrels of con- 
densate would be produced which is 38.6 per cent of the liquid 
originally in the reservoir. Should cycling while continuing 
to market gas be adopted, 26,876,000 million barrels of con- 
densate would be recovered. This is possible because increased 
condensate recovery while the reservoir pressure is still rela- 
tively high is made possible by the larger quantity of gas 
passing through the separators. Under the pressure maintenance 


program nearly all of the condensate, 38,090,000 barrels, can 
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Cumulative Condensate Produced - Millions STB 


Figure 9 
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be expected to be recovered, since the retrograde condensation 
caused by drop in reservoir pressure will not occur. 

A summary of the economic comparison of the three alter- 
native studies is included as Table I. From this summary it 
can be seen that the normal cycling program has the greatest 
present worth to an investor. This results from maximizing 
the value of gas in place by making use of the available mar- 
ket and combining this income with increased condensate recov- 
ery through retrieval at elevated pressure. Earlier retrieval 
of this condensate further contributes to its present net worth. 

Pressure maintnenace results in the maximum ultimate 
recovery of both the residue gas and the condensate. The post- 
ponement of gas sales until the economic limit of the cycling 
of the injection gas is reached, however, makes this procedure 
unattractive economically. Only if there was no market for 
gas at the present time, or if the available market could be 
satisfied by gas from another pool, could investment in pres- 
Sure maintenance prove feasible. In this event, pressure 
maintenance would be the logical choice since the gas would 
have no present value for any of the alternatives and extra- 
neous gas injection to maintain pressure would result in 
valuable recovery of condensate not otherwise obtainable. 

The pressure depletion of the reservoir proves to be the 
least attractive program both from the conservationist viewpoint 


of quantity of hydrocarbons recovered and from the investor's 
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TABLE I 


ECONOMIC SUMMARY 


Pressure Pressure Normal 
Depletion Maintenance Cycling 


Project Life 50 years 80 years 50 years 
Condensate Production, 

M-Barrels 15,430 38,030 26 676 
Total Sales Gas Production, 

MMMSCF 492.7 492.7 492.7 
Ultimate Condensate 

Recovery, % of original 35.6 975710 2 
Gross Working Interest 

Income, M-$ 114,934 17255303 142,860 
Operating Expense, M-$ 4,100 16 367 11,760 
Operating Income, M-$ 110,834 159,000 1 3 0 
Capital Investment, М-5 --- 4,000 5,000 
Total Cash Flow, M-$ 110,834 155,000 126,100 


Present Value, M-$ 38,750 30,000 47,784 
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approach to maximizing profit from the invested dollar. From 
market data available pertaining to recovery from this field, 
some steps toward increasing condensate recovery by gas cycling 
should be taken at the earliest possible time. Since the res- 
ervoir is entering the retrograde phase now, time is of the 
essence if losses from this phenomenon are to be prevented. 

Graphs of the Annual and Cumulative Cash Plans versus 
time are shown for the three alternatives studied in Figures 10, 
mee and 12. 

Capital expenditures required for gas cycling are given 


in Table VIII, page 72. 


Cumulative Cash Flow - Million $ 
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Figure 10 


Annual and Cumulative Cash Flow 
Straight Pressure Depletion 


at Current Rate 
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Annual Cash Flow - Million $ 





Years 





Cumulative Cash Flow - Million 


Figure 11 


Annual and Cumulative Cash Flow 


Pressure Maintenance 
with 


Gas Blowdown 





Annual Cash Flow - Million $ 
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Figure 12 


Annual and Cumulative Casu Flow 


Normal Gas Cycling with 


Pressure Depletion 
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CHAPTER VI 


RECOMMENDATIONS 


Based on the conclusions determined by the investigation 
of the supplementary condensate recovery programs studied, it 
is recommended that the following course of action be under- 
taken by the field operators: 

l. Organize an active operators committee for the 
express purpose of working toward the unitization of the 
reservoir. 

2. Charge the technical subcommittee with the respon- 
sibility of a detailed engineering design for a successful 
normal cycling program. 

3. Charge the appropriate subcommittee with the prepa- 
ration of unitization and operating agreements. 

4, Dedicate the manpower necessary to complete the 
engineering design at the earliest possible time to gain the 
advantage of the highest reservoir pressure attainable during 
recovery. 

5. Commence normal cycling of the residue gas to in- 


crease condensate recovery at the earliest practicable date. 
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Qsc 


DEFINITION OF TERMS 


Water influx 


Effective water compressibility (taking into account 
rock expansion into pore space because of geostatic 
compression of the reservoir rock 


Equivalent field radius 

Cumulative wet gas produced (SCF) 

Initial gas in place (SCF) 

Permeability in darcies 

Porosity 

Viscosity of water 

Effective pay thickness 

van Everdingen and Hurst water influx constant 
Flow rate for well tests | 
Temperature of formation (degrees Rankine) 
Viscosity of gas 

Slope of drawdown curve (psi/cycle) 

Reservoir pressure (psia) 

Gas formation volume factor 

Gas deviation or compressibility factor 
Subscript denoting initial value 

Reservoir pressure decrement 


Infinite aquifer values of dimensionless water influx 
for values of dimensionless time, tp 


Dimensionless time 
Water influx constant 
Real time 


Bulk volume of the reservoir 
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Planimeter 


Area 


(sq. in. 
Ao 21.60 
А, 15. 15 
A, 0265 
А, 19.23 
А, 8,44 
As 6.80 
AL 5. 95 
А. 552 
Аз . 14 
Ag 0 

1 e - 


TABLE II 
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CALCULATION OF INITIAL GAS IN PLACE BY 
VOLUMETRIC METHOD 


) 


Area 
(acres) 


1955 


46 


1391.18 


1161 


32 
1136 
624. 
491. 


47 
12 


0 


61 
39 
02 
42 
27 
75 
86 


Ratio 
or 
Areas 


20 
59 
zo 
-83 
‚80 
. 79 
. 096 
227 
200 


91.827 acres 


Pyramidal formula: 


Trapezoidal formula: 


Sw 


5g 


= 
— 


02 ۷ ٥٣٢ 


ДУ в = 


Interval 


J 


20 
20 
20 
20 
20 
20 
20 
20 
15 


AV, 


7758 x Vg x ó x (1- Sy) x Be 


7758(.28)(.8)(1862.2)(127,825) 


„286 


20 


.000537 bbl/cu.ft. 


Calculation of Flow Capacity 


Khe = 


(m)(P) 


(0.819) (9) (т) (2) (из) 


ог 1862.2 


ДУ 

: 7ء‎ acre-ft. 
Trap 33746 
Trap 25525 
Trap 21010 
Trap 17144 
Ттар 13994 
Trap 111357 
Руши 4613 
Pyr. 569 
Ру 6& 

h -= 

а еи Í 

h 
ü (0۳ Я ee 


O Sos 


SCF/bbl 








TABLE II--Continued 


_ (0.819)(9400)(700)(1.327)(.034) Е 


(34.5) (8000) 


= ,050 darcies 


. 580 3 вя 


60 


С 





1/61 
2761 
3/61 
4/61 
5/61 
6/61 
7/61 
8/61 
9/61 
10/61 
11/61 
12/61 


1961 


1/62 
2/62 
3/62 
4/62 
5/62 
6/62 
7/62 
8/62 
9/62 
10/62 
11/62 
12/62 


1962 


1/63 
2/63 
3/63 
4/63 
5/63 
6/63 
7/63 
8/63 
9/63 
10/63 
11/63 
12/63 


1963 


TOTAL 


MONTHLY PAST PRODUCTION 


MCF GAS 


76,227 
250,783 
11 
409,576 
538, 319 
543, 289 
388, 364 
263,335 
217,357 
433,201 
553,913 
729,040 


4,614,638 


750,948 
696,327 
758,887 
ТАШ ٦ 
756, 234 
618, 763 
675,888 
696, 371 
721,975 
712 
787,513 
788,886 


8,705,775 


723,591 
626,203 
697,915 
610,953 
514,623 
358,040 
625,669 
620,714 
665,982 
666,848 
686,364 
691,439 


7,470,341 


20,790,754 


TABLE III 


BBLS DISTILLATE 


3, 733 
15,098 
16,196 
27,841 
42,100 
42,403 
28,872 
18,530 
15,133 
30,583 
39,992 
54,549 


335,030 


54,465 
52,138 
57,031 
56,523 
55,536 
45,731 
52,441 
47,460 
49,439 
49,682 
60,470 
60,725 


641,641 


61,729 
49,738 
51,312 
49,704 
40,434 
31,177 
47,156 
49,127 
45,531 
46,491 
43,249 
49,899 


565,547 


1,542,218 


61 


BBLS WATER 


8 
368 
399 
130 

86 
282 
143 

35 

84 
138 
300 
217 


2,190 


279 
112 
&65 
&20 
434 
322 
2139 
465 
348 
496 
420 
527 


4,607 





TABLE IV 


PRESSURE DRAWDOWN-BUILDUP TEST 


Offshore, Louisiana 


Pressure Drawdown 


Elapsed 
Date Time n bot сока Remarks 
Press. Psig. 
Days 
3-24-64 9:00 A.M. Wound 72 hour clock 
2 Engaged clock & stylus 
9:15 Pressured lubricator 
9:35 8236 Bomb placed at 11,745' 
10:45 0 8236 Started opening well 
10:52-1/2 .0052 7943 Well flowing 9400 MCF 
11:00 .0104 1292 
11:07-1/2 .0156 7884 
LIND . 0208 7884 
11:30 200985 7876 
11:45 А 7876 
]2:00 noon 0521 7867 
КОЧО Р.М. „0729 7867 
1:00 .0938 7859 
15730 .1146 7867 
2:00 254 7867 
2:30 1507 7867 
3:00 17 7867 
3750 „1979 7867 
4:00 22137 7867 
4: 30 . 2396 7867 
5:00 · 2604 7867 
530 22512 7867 
6:00 3021 7867 
6:30 3229 1099 
7:00 13437 7859 
0ا7‎ . 3646 7859 
8:00 ‚854 2 
S30 .4062 7859 
9:00 +01 7859 
9: 50 . 4479 7859 
10:00 . 4687 7851 
1030 .4896 7851 
11:00 EA Zoo 
2 0 22 7851 


12:00 *midnt.. 25521 7842 
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noon 


TABLE IV--Continued 


Pressure Drawdown 


B Bottom Hole 
m. Press. Psig. 
Days 
20029 7842 
„5937 7843 
.6146 70925 
25957 7025 
. 0962 2917 
2607 | 7613 
¿0979 7809 
. 7187 7809 
. 7396 7800 
. 7604 7800. 
mo 7800 
. 8021 7192 
. O 229 92 
2940 37 7292 
.8646 7188 
. 8854 7788 
. 9062 7188 
‚9271 778& 
. 9479 7784 
. 9687 7784 
29390 7784 

1.9100 7180 

1.0310 7780 

120520 7780 

150730 7780 

1.0940 7776 

1.1150 7776 

21250 7776 

ھ72 1.1560 

1.1770 TZ 

1.1980 7/67 

1.2190 7167 

1.2400 7767 

1. 2600 7763 

122510 / 759 

1. 3020 IDO 

173230 7735 

1.3440 1199 

1.3650 7191 

1. 3559 / 751 

1.4060 7746 


Remarks 


о 





TABLE IV--Continued 


Pressure Drawdown 


Elapsed 
Date Time Time poten ее Remarks 
Press. Psig. 


Days 


3-25-64 9:00 P.M. 174270 7746 
(Cont.) 9:30 1.4480 7742 
10:00 1.4690 7742 
10:30 1.4900 7738 
11:00 1.5100 1156 
230 1.5310 7734 
E2700 midnt. 1.5526 7734 
3-26-64 12:30 A.M. 1.5730 7734 
1:00 1.5940 7734 
230 1,6150 7734 
2:00 1.60 7734 
0۵ھ‎ 1.6560 7730 
BD 1.6770 7730 
3:30 1.6980 7125 
4:00 1.7190 7725 
290 1.7400 7725 
5:00 1.7600 7721 
5:30 1. 7810 142] 
0:00 1.8020 ТТІ? 
6 : 30 1.8230 71 
7.30 1.8650 2717 
8:30 1. 9060 7709 
9: 30 1. 9480 7709 
10:30 1. 9900 7700 
11:30 2.0310 7692 
12:30 P.M. 2.0520 7688 
1:30 2.0730 7683 
2230 2.1560 7680 
3:00 21770 7680 
4:00 2.2190 10795 
5:00 2.2600 76075 
6:00 2.3020 7671 
7:00 2.3440 7671 
8:00 2,3850 7671 
9:00 24270 7671 60 hrs. Shut well in. 





TABLE IV--Continued 


Pressure Buildup 


Elapsed Time 


Surface DWT 


Bottom Hole 


65 


Shut well in. 


Date Time Hours Days Press. Psig. Press. Psie ٣٦٣ 

3-26-64 9:00 P.M 0 0 4745 7671 
9:01 - .0007 6480 -- 
9:02 - .0014 6510 -- 
9:03 - .0021 6530 -- 
9:04 - .0028 6535 -- 
9:05 - 1003834 6542 -- 
9:07-1/2 125 20052 6547 7960 
9:10 167 .0069 6551 -- 
9:15 250 .0104 6552 7977 
9:22-1/2 370 .0156 6556 6985 
9:30 . 500 .0208 6555 7989 
9:37-1/2 620 . 0260 6555 7994 
9:45 120 ‚0312 6554 7994 
10:00 1.000 . 7 6554 7994 
10:15 1,250 0321 6552 7994 
10:30 1.500 ‚0625 6551 799& 
10:45 15557510; .0729 6550 7994 
11:00 2.000 Об 6548 7994 
11:30 2.500 . 1042 6548 7994 
12:00 A.M. 3.000 . 1250 6548 7994 
17:50 3500 . 1458 6548 7994 
1:00 4.000 . 1667 6549 7994 
1:30 4.500 25375 6550 7994 
2:00 5.000 -2083 6552 7994 
230 5.500 ‚2292 6555 وھ‎ ۷ 
3:00 8.000 ‚ 2400 6562 7985 
3:30 6.500 ‚2708 6566 1977 
4:00 7.000 .2917 6570 2977 
4: 30 7.500 23125 6575 7977 
5:00 8.000 223333 6578 7977 
3:30 8.500 ngon 6579 1973 
6:00 9.000 . 3750 6580 7968 
6:30 9.500 . 3958 6583 7968 
7:00 10.000 .4167 6585 7968 
7230 10.500 53175 6586 7968 
8:00 11.000 245089 6586 7968 
8:30 11500 . 4792 6586 7968 
9:00 12.000 0.5000 6586 7968 


Started out 
or hole. 
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TABLE IV--Continued 


Static Pressures 
15,000# Element 


Bottom Hole 


Depth Press. Psip. AP Gradient (Psi/ft.) 
0 6627 -- = 

10,500 7993 ще = 

11,000 8018 25 2050 

11,500 8068 50 .100 

11, 780 8102 34 ‚121 


11,835 








Year 


Q «000-109 ی‎ о м н 


Condensate 


Bbls 


617490 
585247 
547395 
Б 7187 
482908 
456178 
429448 
409448 
389447 
272718 


356082 
342718 
333745 
218231 
311502 
300 380 
295987 
289259 
282623 
229137 


273744 
271501 
267015 
262623 
260380 
258136 
255893 
205993 
253650 
251407 


STRAIGHT DEPLETION OF PRESSURE 


Wet Gas 
MMSCF 


9346 
9346 
9346 
9346 
9346 
9346 
9346 
9346 
9346 
9346 


9346 
9346 
9346 
9346 
9346 
9346 
9346 
9346 
9346 
9346 


9346 
9346 
9346 
9346 
9346 
9346 
9346 
9346 
9346 
9346 


TABLE V 


PRESENT WORTH 


AT CURRENT RATE 


Gross 
Income 


3052287 
2973027 
2881280 
2794409 
2625995 
2658742 
239 202 
2543778 
2495943 
2455130 


2414540 
2381940 
2360050 
2322200 
2305784 
2278650 
2267935 
2241420 
2235230 
2224385 


2213668 
2208195 
2137259 
2155599 
2181064 
2175589 
2170117 
2170117 
2164646 
2199172 


Oper. 
Costs 


124439 
120005 
114801 
109904 
105934 
192259 
28599 
95834 
93084 
90784 


88494 
86658 
85424 
83291 
82364 
80836 
80234 
79307 
78394 
71716 


77174 
76864 
76249 
75644 
75334 
75028 
74719 
74719 
74411 
74102 


Undeferred 
Net Income 


2927848 
2853622 
2766479 
2684505 
2618061 
2556484 
2494949 
2447944 
2402859 
2364346 


2326046 
2295282 
2274626 
2238909 
2223420 
2197814 
2167099 
2172213 
2156956 
2146607 


2136494 
ОТЕ 
2121001 
2110889 
2104730 
2100561 
2094398 
2025598 
2090235 
2085070 


Deferred 
Net Profit 


2842940 
2616770 
2390240 
2190556 
2013290 
1856000 
1790040 
1581979 
1465640 
1359509 


1263040 
1175185 
1098645 
1020940 
953850 
890150 
835700 
784170 
733360 
689060 


647360 
609460 
572670 
563165 
3205375 
&7&725 
448415 
423875 
397145 
375310 
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Year 


TOTAL 


Condensate 


Bbls 


250380 
249258 
248136 
244770 
244770 
244770 
244770 
244770 
244770 
244770 


244770 
244770 
244770 
244770 
244770 
244770 
244770 
244770 
244770 
244770 


15429856 


Wet Gas 
MMSCF 


9346 
9346 
9346 
9346 
9346 
9346 
9346 
9346 
9346 
9346 


9346 
9346 
9346 
9346 
9346 
9346 
9346 
9346 
9346 
9346 


467300 


TABLE V--Continued 


Gross 
Income 


2146667 
2153939 
21511 0 
2142980 
2142980 
2142980 
2142980 
2142980 
2142980 
2142980 


2142980 
2142980 
2142980 
2142980 
2142980 
2142980 
2142980 
2142980 
2142980 
2142982 


Oper. 
Costs 


73961 
73807 
13552 
73190 
72190 
73190 
73190 
73100 
70 
73190 


73190 
73190 
73190 
73190 
7/0 
70 
73130 
73190 
7/0 
73190 


Undeferred 
Net Income 


2082706 
2080123 
2077557 
2069790 
2069790 
2069790 
2069790 
2069790 
2009720 
2069790 


2069790 
2069790 
2069790 
2069790 
2069790 
2069790 
2065790 
2069790 
2069790 
2069790 
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Deferred 
Net Profit 


347810 
330740 
312710 
295910 
29595] 
262860 
246300 
29 9990) 
219400 
206979 


195180 
184830 
173860 
164340 
154820 
146334 
140124 
130600 
123150 
115910 


38749249 











TABLE VI 


PRESENT WORTH 
PRESSURE MAINTENANCE WITH GAS CYCLING 
40 MMSCF PER DAY 
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Condensate Wet Gas Gross Oper. Undeferred Deferred‏ مات 
Bbls MMSCF Income Costs Net Income Net Profit‏ 
)2059800 -( 1998150 441500 2439650 1000000 1 
)227497-( 1998150 441500 2439650 1000000 2 
1498900 1998150 441500 2439650 1000000 > 
1630490 1998150 441500 2439650 1000000 4 
1536600 1998150 441500 2439650 1000000 5 
T50 1450660‏ 19983 441500 2439650 1000000 6 
D 1000000 2439650 441500 1998150 1368700‏ 
1290800 1998150 441500 2439650 1000000 8 
1218900 1998150 441500 2439650 1000000 9 
1149000 1998150 441500 2439650 1000000 10 
1085000 79938350 441500 2439650 1000000 11 
T50 1023000‏ 1998 441500 2439650 1000000 12 
965100 1998150 441500 2439650 1000000 13 
911160 1998150 441500 2439650 1000000 14 
857200 1998150 441500 2439650 1000000 15 
809250 1995150 441500 2439650 1000000 16 
7103500 1998150 441500 2439650 1000000 17 
721300 1998150 441500 2439650 1000000 18 
е) 1000000 2439650 441500 1998150 679370‏ 
641400 1814000 430500 2244500 920000 20 
605440 1744925 426375 2171300 890000 21 
571500 1675650 422250 2098100 860000 22 
259500 1606775 418125 2024900 830000 25 
507500 1560725 415375 1976100 810000 24 
479560 1491650 21.1250 1902900 780000 25 
451580 1422575 407125 1829700 750000 26 
427600 1376525 404375 1780900 730000 27 
403630 199055 401625 1732200 710000 25 
379650 1284525 398875 1683400 690000 29 
359670 1238475 396125 1634600 670000 30 
UBTOTAL 27640000 24325760‏ 
RESERVOIR BLOWDOWN‏ 
5687293 102179117 3446150 105025267 482.3 10450000 80 
TOTAL 38090000 482.3 3001505‏ 


, 
Р! 





Near 


О р о м O سم یم نن ےہ ی‎ 


| 


Condensate 


Bbls 


1605475 
1521642 
1423230 
1330646 
1255560 
1186062 
1116565 
1064565 
1012562 

969067 


925813 
891067 
867737 
744660 
688420 
632300 
588720 
545253 
510700 
477280 


273744 
271500 
267015 
202673 
260 380 
259106 
255893 
255893 
253650 
241407 


Wet Gas 
MMSCF 


9345. 
9345. 
9345. 
9345. 
9345. 
9345. 
9345. 
9345. 
9345. 
9345. 


On ON O) ДУ О ON ON ON ON ON 


9345. 
9345. 
9345. 
9345. 
9345. 
9345. 
9345. 
9345. 
9345. 
9345. 


ON O9 N ON ON TS ON 


9345. 
9345. 
9345. 
9345. 
9345. 
9345. 
9345. 
9345. 
9345. 
9345. 


Ov, OV OV ON Cn Ov CON OY ON ON 


TABLE VII 


PRESENT WORTH 
GAS CYCLING WITH PRESSURE DEPLETION 


Gross 
Income 


5462600 
5258100 
5018000 
4792100 
4608950 
4439400 
4269860 
4143000 
4016130 
3910000 


3804500 
37 19700 
3662800 
3362540 
32 25330 
3088420 
2982100 
2876060 
21700 
2710244 


2215070 
2208200 
2197250 
2186540 
2181060 
2175600 
21701 Ó 
2170129 
2164650 
2149170 


GAS SALES AT CURRENT RATE 


Oper. 
Costs 


564285 
552760 
539290 
526500 
316175 
506620 
&97065 
399015 
482765 
476785 


470835 
466055 
462845 
445925 
439193 
430475 
424485 
418535 
413755 
409165 


ТҮЙІП» 
76865 
76265 
75645 
7ھ‎ 
73029 
74720 
74720 
74410 
74105 


Undeferred 
Net Income 


4898315 
4705340 
4478770 
4265600 
4092775 
3952780 
900259» 
3653086 
5555265 
3433215 


3333665 
3253645 
3199955 
2916615 
2787137 
2657945 
2557615 
2456525 
2378005 
2301079 


2136495 
1231333 
2120985 
2110895 
2105725 
۲ 2ھ" 
2095400 
2095400 
240 2090 
2085065 


Deferred 
Net Profit 


(-243735) 


4071065 
3869660 
3480730 
3147350 
2855200 
2584360 
2359909 
2155559 
1974100 


1810180 
1665870 
1545600 
1330000 
1195709 
1076467 
977000 
887170 
808521 
738810 


647360 
609560 
572670 
536170 
3205575 
474730 
448415 
2523270 
397145 
95910 
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TABLE VII--Continued 


Year Condensate Wet Gas Gross Oper. Undeferred Deferred 
Bbls MMSCF Income Costs Net Income Net Profit 

| 250380 9345.6 2156070 7/3965 2082705 347810 
12 249248 9345.6 2153930 73805 2080125 330740 
35 248136 9345.6 2151190 73665 2077555 3200 
34 244770 9345.6 2143000 73190 2069810 ZI 
35 244770 74926 2143000 (5190 2069810 277550 
36 244770 9345.6 2143000 73190 2069810 262670 
37 244770 9345.6 2143000 73190 2069810 246 300 
38 244770 9345.6 2143000 73190 2069810 235900 
39 244770 9345.6 2143000 ل7‎ 2069810 219400 
40 244770 9345.6 2143000 73190 2069810 207000 
4l 244770 9345.6 2143000 73190 2069810 195180 
42 244770 9345.6 2143000 751090 2069810 184800 
43 244770 9345.6 2143000 73190 2069810 173860 
44 244770 9345.6 2143000 990 2069810 164340 
45 244770 9345.06 2143000 ono 2069810 154820 
46 244770 9345.6 2143000 73190 2069810 146340 
47 244770 9345.6 2143000 73190 2069810 140130 
48 244770 9345.6 2143000 73190 2069810 130600 
49 244770 9345.6 2143000 73190 2069810 122150 
50 244770 9345.6 2143000 73190 2069810 115900 
TOTAL 26876429 467280 47784138 





TABLE VIII 


CAPITAL EXPENDITURES 


Pressure Maintenance - 


Recomplete 9 
New Well 1 
Recomplete 0 
New Well 5 
Recomplete 4 
Well 6 


platform 


Fabrieagion 
Installation 


Compressors 


Fabricat ron 
Installation 


Pipe 


High pressure line 


Sales gas line 


Separation Facilities 


Cycling without gas sales 


5 70,000 
600,000 
70,000 
600,000 
70,000 

0 


320,000 


80,000 


1,116,000 


99 0 


25,000 
100,000 


Eds tional Plant Required for Normal Cycling 


New well 


Sales gas line 


5577902000 
300,000 


Total for normal gas cycling 


$1,410,000 


$ 400,000 


$1,206,000 


$ 125,000 


9 971,500 


$4,000,000 


$1,000,000 


$5,000,000 


Да 





APPENDIX I 


MATERIAL BALANCE CALCULATION 


E 
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GAS RESERVE CALCULATION PROGRAM 


Dimension 
Common 





Read Input: Data for 
Subroutine Z calc constants 


Read Input Data Lor ` 
Q funct Table 


Read Input: Data for 
Reservoir Properties 


Read Input: Gas, I 
Distillate, and Water Produced 


Print Out Reservoir Properties 


Read in Pressure Decrements 


Initialize Constants 
XW = 0 YW = 0 


Xu YY : 
WW = 0 


ú 
4 





C = .00515 (Temp) 


с _ 1858065 GO (T) 
2 (а) (@) (т) (т) 


6 = Р(1) 


Calculate 7 





| - 
ET = C x 278 | 
у= 
ier (0) 
X = Ü ' 


— — — — Dé2, K=2, Inc, l 
L = k-1 

A P(K) = Р 1 - P. 

Calculate Influx WD 


Calculate Z 
С х 7/5 


WI k) = WD 
|W = WD | 
Š = P 

Be = 





~ Sow (k) 


5 


XY (С 09) C) * XY 
YW = (B -B ;) (Wp) + YW 
WW = (W) Gr) + XX 

и = 


у = he 
(G М 
g D 
- у^ + үү 
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Calculate Gas in Place: 


с = QU)QW) - QW) Qm) 


— (WW)(YY) - (YW) (YW) 


Calculate Water Influx Constant: 


XW - (G)(YW) 


AA = WW 





PRINT G, AA 
WE = 0 
eo 


Print Initial Conditions 
Time, Pressure, AP, Influx 


DOCG EF S 


LK 


Calculate Water Encroachment 


WE = (AA) (WI (o? 


Calculate Gas 
- WE 














сБ 
ОЭ, 
А B8 (1) 






Calculate Deviation 


PCent = рана 






= ے‎ | Continue 








800 


101 


2ء 


1000 


201 
ص9 


202 
1020 


203 


204 
1050 


102 


300 
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PROGRAM GASRES 

DIMENSION Al6s6) 9B(6) 51 (381) 501381) sDELP (20) sP120)5GP120)sDPiZ0 
DIMENSION WP(20) sGPW(20) »М1 (20) „ОЕГВО (20) 
СОММОМ А,В, ТЕМР „РС „ТС ,Т „СО »К „ОЕ! Р 

PRINT 800 

FORMAT(1H1) 

PRINT 101 

BSEMAT(44H RESERVOIR COMPUTER PROJECI---—-ENERGY BALANCES 
DATA INPUT FOR SUBROUTINE ZCALC CONSTANTS 
DO 1000 J=1>691 

DO 1000 I=1>691 

READ 200sA(I J) 

FORMAT (13X9F14.11) 

CONTINUE 

ОО 1010 1=1,6,1 

READ 201,B(T) 

FORMAT(F14.11) 

CONTINUE 

DATA INPUT FOR SUBROUTINE QFUNCT TABLE 

ОО 1020 141>381»3 
БЕСБФОООЭТТСТТІЭОСТТЭТСТЕТТУОГТЕТОЭТТТЕ TAE) 
ФЕМАТ(Ғ12.2>Ғ12.3%Ғ12.25Е12.3»Ғ12.2»Ғ12%323) 
CONTINUE 

BERTA INPUT FOR RESERVOIR PROPERTIES 

READ 203:۸ 1۷۴۱ 

READ  )۷ 

READ 203,RI 

READ 203»ПЕРТН 

READ 203, TEMP 

READ 203:56] 

READ 20350 

READ 203%9AH 

READ 203956 

READ 2039U 

READ 203»AK 

BEAD 2039TIME 

READ 203.568 

READ 203sSWC 

READ 203•Р1 

КЕРІ) 209ЭРС 

READ 203,ТС 

КЕАО 203, Ви 

READ 2039CPI 

FORMAT(F16.7) 

ІМС-СРІ 

DATA INPUT FOR RESERVOIR PRODUCTION 

DO 1030 T=1sINCol 

READ 204>GP(I)»DP(I)>WP(T) 

FORMAT (3F20.0) 

CONTINUE 

Pee) 102 

FORMAT(19H RESERVOIR DATA FOR3/) 

TS FOIS UOT TOR RESERVOIR PROPERTIES 
PRINT 300sINC 

FORMAT(45H INC sCURRENT NUMBER OF PRODUCTION 1ТМСКЕМЕМТ 5=» 14) 
PRINT 301 sATIWEL 





E 
502 
203 
304 
305 
306 
307 
EOS 
309 


210 


314 
215 
316 
ЗИТ 
802 

1 

205 
1040 

a 


103 
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FORMAT(48H MAX. NUMBER OF PRODUCING WELLS AT ANY ONE ТІМЕ-»Ғ6.1) 
PRINT 302 5BW 

FORMAT(50H BWsFORMATION VOLUME FACTOR FOR WATER IN BBLS/STB=>F6+3) 
PRINT 303,CONV 

БОБИЛ ЗАН СОМУ, 5СҒ OF VAPOR PER BARREL OF DISTILLATE ST EN) 
PRINT 304,RI 

ESEMATUGBH RISRADIUS TO INITIAL EDGE OF RESERVOIR IN TEES 
ШОШЫП 305>5DEPTH 

۰... ۱9م‎ DEPTHsDEPTH OF RESERVOIR IN РЕЕТ-„ЕВ.! ) 

PRINT 306 و‎ ٘ ۳۳ 

SORMAT(SOH TEMPS TEMPERATURE OF RESERVOIR IN DEGREES RANE TE Sy 
PRINT 307 9PC 

menial’ 30H PC»sPSEUDOCRITICAL PRESSURE IN PSIA=,F6.1) 

ӘКІМІ 308 „ТС 

SI JAT(S5SOH TC |+ PSEUDOCRITICAL TEMPERATURE IN DEGREES RANK TNE SSF orii 
НАТ 309sSGI 

FORMAT(40H SGI» INITIAL GAS SATURATION AS FRACTION=sF5.3} 

ВВГ 0 

FORMAT(24H OsPOROSITY А5 ЕКАСТІОМ-»Ғ5.3) 

РКІМТ 311»АН 

FORMAT(43H AHsAVERAGE THICKNESS OF RESERVOIR IN FEET=oF6e1) 
INT 312sSG 

monn! (44H SGsSPECIFIC GRAVITY OF GAS RELATIVE TO ATR STI 
ВЕТ 313.0 

ШЫ МАТ(Л6Н UsVISCOSITY OF RESERVOIR WATER IN CENTIPOTSES Hon 
PRINT 314,AK 

FORMAT(47H AKS.WATER PERMEABILITY OF RESERVOIR IN DARCIESz,F7.e4) 
РЕТАТ 315,TIME 

FORMAT(33H TIMESTIME INTERVAL USED IN DAYS=F7.2) 

ERTNT 316>SGR 

FORMAT( 41H SGRsRESIDUAL GAS SATURATION AS FRACTION=5F5-3) 

ee? 317>SWC 

FORMAT(42H SWCsCONNATE WATER SATURATION AS FRACTION=5F 503} 

PRINT 802 

FORMAT( 1H 9/7) 

fee INPUT FOR RESERVOIR PRESSURE 

DO 1040 I=1>INC>1 

MEAD 2059P(1) 

ШЕБГАТ(Р10.3) 

CONTINUE 

MAIN BODY OF PROGRAM WITH CALCULATIONS 

INT 103 

FORMAT (16H TIME PRESSURE » 9 ۸۰۹۰1۳12 10 و۸۴۰‎ 100۳166190۷۳۰۸ 5 6 ٣٣٣٢ БС ٥ 


1s14HRES. WITH.sBBL) 


DELP(1)=0.0 
XW=0.0 

YW=0.0 

XY20,0 

YY=0.0 

WW=0.0 
С-0.00514590%ТЕМР 
DELTD21858064,.0X*AK*TIME/(U*O*RI*RI) 
S=P (1) 

CALL ZCALCI(S>2) 
BGI=C*Z/S 

7-04 0 
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L=0 
Х-0.0 
PRINT 4009L9S3sZ5BGI 9Y 9X 
ESOCSOMEIAT(IISS2XSFIO0Q3$.2X*sET4A8s2X9*ET4.8.32XSEI4. 8B$2X VETA4S B) 
DO 2 K=2sINCs1 
L=k-1 
DELP(K)zP(K-1)-P(K) 
WALL INFLUX(DELTD sWD ) 
WI(K)=WD 
W=WD 
S=P(K) 
CALL ZCALC(SsZ) 
BG=C*#Z/S 
GPW(K)=GP(K)*BG+WP(K) *BW+DP(K) *CONV*¥BG 
X=GPW(K) 
Y-BG-BGI 
DELBG(K)=Y 
XW=X*W+XW 
ХУ-ХХУ+ХУ 
УИЗУКИ+ УИ 
mer =¥ #Y +Y Y 
WW < ل۶ إ۷‎ +۷۷۱۷۷ 
PRINT 400sL9S9ZsBGsY 9X 
EECONTINUE 
G-sCOXY*WW)-OXWYWy )/ UCWw*YY)—CYW*YW) ) 
AAZUXW-G*YW) / WW 
PRINT 401»С»АА 
401 FORMAT(1H »//s»14H L.M.S. G»SCF=sE14.8»/»23H WATER INFLUX CONSTANT= 
29E14e83//) 
PRINT 104 
104 FORMAT(16H TIME PRESSURE s6Xs5HDEL Ps3Xs14HWATER INFLsBBLs2Xs14HP 
КЕКТООТС G CHK »2Xs12HPERCENT DEV.) 
МЕ-0.0 
(-0 
PRINT 402s5LsPIsDELP(1) ۶۳۴۳ 
402 FORMATLI5»2X»F10.3s2X»F10.3»2XsE14.8) 
DO 3 К= 2» ІМС» 1 
L=K-1 
WE=AA*WI(K) 
GC=(GPW(K)-WE)/DELBG(K) 
PCENT=(GC=6)*100.0/6G 
РКТМТ 403,| »Р (К) „ОЕ:Р (К) + МЕ» СС» РСЕМТ 
4803 ҒОКМАТ(1І5»2Х>Е10.3»>2Х>Е10.23»2Х>Е14.8»2Х>Е14.8%2Х»ҒЕ10.4) 
3 CONTINUE 
PRINT 801 
801 FORMAT(1H s/////) 
PRINT 800 
CONTINUE 
Bo TO 1 
END 
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SPERSUTTNE ZCAECUSSZ) 

DIMAS ۱۹3۲ Al6s6)»B(6)>10381)501381)>DELP (20) sP 2053562207, 027203 
DIMENSION PX(6)sPY(6) 

COMMON AsBsTEMPsPCsTCsT sQsKsDELP 

А ОХТЕМРИТС) -4.0}/1.9 

posu sagt) 

2(2) ۶۰۲۷ 

moo) =8(3)*((3.0*Y*Y)=—-1.0) 

Ey) =B(4)4#((5-0*Y#Y) - 39 ۷پ‎ 

fee )=68(5)*(((35.0*Y*Y )—-30,0)*Y+7+3.0) 
beeen) =816)*((((63.0*Y*Y )-70.0) КУХУ) + 15 ,„ОРкУу 
PX(1)=B(1) 

PX(2)=B(2)*xX 

Meer) =81(3)*((3.0*X#X)-1.0) 
PX(4)=B(4)¥((5.0*X*¥X)-3.0) *X 
Px45)=B(5)*((((35,08X#X )-30.0) *X*X | +3. 0) 
PX(6)=B8(6)*((((63.0%X*X )~70.0) XX*X)4-154,0) *X 
APP=0.0 

DO 1 M=1>691 

DO 1 N=1>691 

APP=(A(NsM)*PX(M)*#PY(N) )+APP 

Z=APP 

RETURN 

END 
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SUBROUTINE QFUNCT(DELT sQT) 
DIMENSION A(626) »B(6) 9T( 381) »>Q(381) sDELP(20) »P (20) » СР(2 0)» рР(20) 
COMMON A»B,TEMPs»PC»TC»T»Q»K»DELP 
TT-DEET 
DO 3 Іг19381,1 
J ([T—T(1)) 1s253 
pee (T —1 )=G[(1—2))/ (T(1-1)=T(/(1=2)) 
al ON CU E A A E AA ٦ 
ПО -2) + (ТТ-Т(1-2)) “001+ ( ТТ-Т 12; | ( ТТ-ТЕТ- 1) 1 2002 
GO TO 4 
QT=Q(1) 
GO TO 4 
CONTINUE 
RETURN 
END 
END GASRES 
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SUBROUTINE INFLUX(DELTDs wD} 
DIMENSION A(656) 9B(6) 9T( 381) 9Q(381) sDELP(20) sP(20) »GP(20) sDP(20) 
COMMON AsBsTEMPsPCsTCsTsQsKsDELP 
WD=0.0 

DO 1 J=1>K>1 

TI=J=1 

DELT=TI*DELTD 

EL GEUNCT(DELT sQT) 

N=K+1=-J ` 

WI=DELP (N) #QT 

WD=WI+WD 

CONTINUE 

RETURN ` 

END 





КЕРЕШ ГЕ COMPUTER PROJECT—-=“ENERGY BALANCE 


NESERVOTR DATA FOR 


INCSsCURRENT NUMBER OF PRODUCTION INCREMENTS- 13 


MAXe NUMBER OF PRODUCING WELLS AT ANY ONE TIME= 2.0 

BWsFORMATION VOLUME FACTOR FOR WATER IN BBLS/STB= 1.030 

٠۰۰٣ع ع0‎ УАРОВ РЕВ BARREL OF DISTILLATE= 768.0 

E [US TO INITIAL EDGE OF RESERVOIR IN FEET= 2190.0 

EEEEHSDEPTH OF RESERVOIR IN FEET= 11450.0 

MEMES TEMPERATURE OF RESERVOIR IN DEGREES RANKINE= 690.0 

APS EUDOCRITICAL PRESSURE IN PSIAS 662,3 

[ieee SEUDOCRITICAL TEMPERATURE IN DEGREES RANKINE= 38765 

SGIs INITIAL GAS SATURATION AS FRACTION= „800 

۰۰۰۰٥٥. ۱۲۷ AS FRACTION= „280 

AHsAVERAGE THICKNESS OF RESERVOIR IN FEET= 64.4 

OS PECIFIC GRAVITY OF GAS RELATIVE TO AIR= .7975 

ШОО COSITY OF RESERVOIR WATER IN CENTIPOISE= 300 

AKsWATER PERMEABILITY OF RESERVOIR IN DARCIES= ,0500 

TIMEsTIME INTERVAL USED IN DAYS= 90.00 

SGR.SRESIDUAL GAS SATURATION AS FRACTION= .300 

SWC>»CONNATE WATER SATURATION AS FRACTION- „150 

TIME PRESSURE А BG»BBL/SCF RES ۷۷۷٣٦ 
0 8838.000 1432981926570 Gre 3. 3 ۵٥ 0 
1 8818.000 1.32803329۴ 068 5.634748162E-04 4407006557202 
2 8765.000 1.32322017Е 00 5.36093049E-04 1,149964894Е 06 
В 8717-000 1.231904095Е 00 5.37281226Е-04 1.04249002Е 06 
4 8657.000 1,31371808E 00 5.38827841E-04 250 273 277200 
5 8582.000 1.30708189Е 00 5.40 765 1ОЗЕ-04 3:0722 TOER 
6 8509.000 1.30063895Е OO 5.42736044E-04 5.12844213E 06 
7 8440.000 1.29456183Е OO 5.44616486Е-04 6.35037990Е 06 
8 8364.000 1.28788050E ОО 5.46728833E-04 17.69974695Е 06 
E 8289.000 1.28129781Е 00 5.48855951Е-04 Зл 1264 ВЗТЕ СБ 
10 6219.000 1.27516200Е ОО 5,508719 758E 02a 7 ЈЕНИ 
0 9155.000 1.26955778E 00 5 ے2 9ء‎ ۲٣29 ۹  ۷'؟+؟؟‎ ЕКТ 
12 8100.000 1.264 14529Е 00 5 54606719٢ - 08۰٦ 

LeMeSe GeSCF=6.06624695E 11 


WATER INFLUX CONSTANT=4.0706801*E 00 


TIME 


بر رم یں می ص0 


فم فم يم 
fN — O‏ 


O SURE 


8818.000 
8764.000 
8717.000 
8657. 000 
8582.000 
8509.000 
8440.000 
8364.000 
8289 «000 
8219.000 
8155.000 
8100.000 


DELTE 
20.000 
54.000 
47.000 
60.000 
5000 
ОО 
69.000 
76.000 
75.000 
OE OE 
64.000 
5.000 


WATER INFLsBBL 


6 707۳ 


1.0311549*E 
4.534 754 7*E 
Яо 
l6 e Oe ais 
2+655 267۴ 
30000022 Ome 
5.2608605*Е 
60090593055 
80 729 7406*E 
1.0 748285 *#E 
1022188892 E 


6. 15121449E 
6,26 332 r 
5.45345143E 
5. 774 76619Е 
5.298680) Е 
6.07865842Е 
6.16642836Е 
6 | 21517592EF 
6. LS 229905 
5.95654287E 
5.98 3 54920Е 
6.11444376Е 


DAE 
LI 
1 
uon 
Га 
11 
T3 
a 
it 
11 
П 
1 





APPENDIX II 


FLUID DATA 





Component 


Carbon Dioxide 
Nitrogen 
Methane 
Ethane 
Propane 
Iso-Butane 
N-Butane 
Iso-Pentane 
N-Pentane 
Hexanes 
Heptanes Plus 


Totals 


HYDROCARBON ANALYSES OF SEPARATOR PRODUCTS 


Separator Liquid 


Mol % 


AND CALCULATED WELL STREAM 


Lig. Vol. % Mol % 


N GQ) = O = = = = Q 
C2 اما‎ 
~ хо 

хо 


со 


| 
© 
© 
O 
© 


Calculated Specific Gravit Air = 1.00 


Eroperties of Heptanes Plus: 


API Gravity 


Specific Gravity 
Molecular Weight 


Basis of Recombination: 


39.9°API @ 60°F 
0.8256 @ 60/60°F 
202 


Мә 
lo 
сос O O O O O F Q Ç O O 


100. 


Separator Gas 


15 


СРМ@ 
.025 PSIA Mol % 
0.67 
0.30 
89299 
3556 
0.349 1. 7 
8 0.43 
0.081 0.35 
0.041 0.17 
0.018 0.13 
02.029 0230 
0.019 2.82 


0.654 100.00 


0.606 


Separator Liquid per MMSCF Separator Gas = 66.07 Bbls 


85 


Well Stream 


GPM@ 


15.025 PSIA 


. 385 
· 144 
SS 
.063 
‚048 
. 1206 

„209 


N O O O O O O 


3.112 


027975 





Relative 
Pressure Volume 
(PSTA) (V/Vsat) 
8680 0.9793 
8600 Original Reservoir 
Pressure - 
8305 0.9978 
8260 Dew Point Press. 1.0000 
8138 1.0070 
7825 1.0255 
2555 1.0442 
6930 1.0910 
6415 1.1382 
5990 1.1856 
5610 172392 
5295 1.2808 
4765 1.3767 
4168 1.5207 
3495 1. 7555 
2675 2. 2801 
2190 2.7208 
1855 3.2067 
1620 3.6864 


PRESSURE-VOLUME RELATIONS OF RESERVOIR FLUID AT 240?F 
(CONSTANT COMPOSITION EXPANSION) 


EFFECT OF GAS-LIQUID 


Dew Point Pressure 
(PSIA) 


8,090 
8,260 
8, 360 


Specific Deviation 
Volume Factor 


(Cu. Ft ШЫ) (2) 


0.04918 0+817 


.05011 | 2782 
-05022 1.2741 
.05057 
„09150 
„05244 
„05479 
56 
.05954 
"06195 
.06432 
06914 
.07657 
.08816 
11250 
"156063 
.16104 
.18513 


> جح ا‎ O Оо со о оо CO осо OO оо соо 


КАТТО UPON DEW POINT PRESSURE АТ 240°Е 


Gas-Liquid Ratio 
(SCF Separator Gas per 
Bbl Separator Liquid) 


20,000 
15,135 
10,000 





| 





“ 
1 


87 


RETROGRADE CONDENSATION DURING GAS DEPLETION AT 240°F 


Reservoir Liquid Retrograde 
Pressure (Bbls/MMSCF of Liquid Volume 
(PSIA) Reservoir Fluid) Percent * 
8,600 Original Res. Pressure 0 0 
8,260 рем Point Pressure 0 0 
7,500 14,09 2.48 
6,500 31.55 52:55 
5, 500 42.33 7.45 
4,500 47.03 0227 
3,000 48.64 8.56 
1,500 &5.57 6.68 


* Percent of Reservoir Hydrocarbon Pore Space 


HYDROCARBON ANALYSES OF GAS AND LIQUID REMAINING IN RESERVOIR 


AFTER DEPLETION TO 1,500 PSIA at 240°F 


Reservoir Gas Reservoir Liquid 
СРМО 
Component Mol % 15.025 PSIA Mol % Liq. Vol. % 
Carbon Dioxide 0.74 
Methane and Nitrogen 92.36 26.74 7.42 
Ethane ۶ 0 3 22 1. 38 
Propane 1.46 0.410 1.74 0.79 
Iso-Butane 0.44 On 147 0.65 0.35 
N-Butane 0531 0.100 0.65 0.34 
Iso-Pentane 0.17 0.063 0.98 0.59 
N-Pentane 0.11 0.041 0.74 0.44 
Hexanes 0.21 0.088 1.99 1.34 
Heptanes Plus 0.50 0.251 63.29 87.39 
Totals 100.00 1.100 100.00 100.00 


Properties of Heptanes Plus: 


API Gravity 2-2 ови ۳۲ 
Specific Gravity 0.8343 
Molecular Weight 125 


ни 





HYDROCARBON ANALYSES OF PRODUCED WELL STREAM (MOL PER CENT) 


Component 


Carbon Dioxide 
Methane and Nitrogen 
Ethane 

Propane 

Iso-Butane 

N-Butane 

Iso-Pentane 
N-Pentane 

Hexanes 

Heptanes Plus 


Totals 


DEPLETION STUDY AT 240°F 


Dew Point 


Fluid 
8260 


90. 20 


Properties of Heptanes Plus 


Specific Gravity 
Molecular Weight 


0.8256 


202 


7500 


0.75 
DIS 
3.41 


196 


Reservoir Pressure - PSIA 


6500 


0.8198 0.8146 


186 


5500 


0.8058 


169 


GPM Content of Produced Well Stream (Gal/MSCF) 


Propane 
Iso-Butane 
N-Butane 
Iso-Pentane 
N-Pentane 
Hexanes 
Heptanes Plus 


Totals 


Deviation Factor "Z" 


Well Stream Produced 


Cumulative Per Cent 


. 385 
. 144 
113 
. 063 
. 048 
. 126 
+0, 


3.112 


сос O O O O O‏ بح 


1.280 


0.349 
. 130 
. 100 


09-315 
0.120 
0.087 
0.060 
0.044 
0.096 
1.249 


1.971 


12.125 


LON 


2322 
LIT 
.084 
. 060 
.044 
.088 
.841 


7 


Ce 0 о> о ооо 





بے[ 


1.054 


poro" 


4500 


051981 


156 


223 
SEI 
.084 
.060 
.041 
‚084 
‚282 


229] 


O O O O O O O 


[= 


02992 


29.84 


3000 


0.7826 


157 


"957 
‚130 
2087 
2056 
.041 
.080 
.368 


nn 


1.119 


تہ ہت © ооо‏ 


0.927 


49.56 


88 


оне и Шо мо uw. 1 ——— s о = 0 вв 


0.7665 


119 


‚410 
.147 
‚100 
‚063 
‚041 
. 088 
251 


‚100 


O OOO 0202 


E 


0992 


13.51 


87i 





89 
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Relative Volume (V/Vsat) 
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PRESSURE-VOLUME RELATIONS OF RESERVOIR FLUID @ 240°F 


Rh 
00 





= 
о 






EJ 


= 
N 


لم 
I‏ 


-- 
© 


2000 4000 
Pressure (PSIA) 





Gallons per Thousand Cubic Feet 


GASOLINE CONTENT OF HYDROCARBONS PRODUCED 


DURING PRESSURE DEPLETION @ 240°F 


Pressure (PSIA) 


92 
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DEVIATION FACTOR "Z" OF GAS PHASE 


DURING DEPLETION @ 240°F 


9 Experimental 


© Calculated 





x 
Ml 
x 
К 

wi 
Lr 
ہے‎ 
A 


со 


ہے 


2 


|. 


8000 
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ü 
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DURING DEPLETION @ 240°F 
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RETROGRADE CONDENSATION DURING 


DEPLETION @ 240°F 
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4000 
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POTENTIAL DISTRIBUTION PROGRAM 
PROGRAM POTMAP 


FAIL = 0.0 







— — — — — pO 29 о 





Conditions 


Set Boundar 
б CONST 7 106000./P. . 
EG 








98 


(-CONST) (Q, .) + (RK, RK 


+ 

elle 
RK + ВК. . 
08) 


С оС ЕЕ" 


і p Ciri j ГЕ) 
T АКК, 


ЕШ + RK. 
1- 


+ 
i jl 1° FÉ jy id 


2 


+ +1) + + 
+ (ВК. "mus Depi 1) + (RK, $ RK .) (P 


i 1,3 j em 


IF 
rror = ABSF (PN-P. .) 


2 
(в) e^ 
| 
| 
| 
| 
| 
| 
Е. 


= 


FAIL = 1.0 
1, 


— — — — -- — — — — — CONTINUE 


PRINT FAIL 
A 
= 0 
| 
| 
| 
Lo — === ===> CONTINUE 


PRINT Output 


Poe 
1, 


END 





40 


42 


ee 


46 


48 


43 


45 


om 


12 


ШЕ 


LS 


PROGRAM POTMAP 
3.۱۶۸510۸ P(33,38),RK(335,38)50133,58) 


00 40 1=1.31 

БЕЛО 625 (P(lsJ).J=1512) 
ВО 42 1=1.31 

о бог • (Р(1»Ј)»2=13,24) 
DO 63 121%+3] 

БЕС) O62s(P(1lsJ)s J=25536) 
DO 46 1=1,31 

READ 41» (Q(I sJ) 9J=1318) 
DO 48 ٢-1 1 

FAD 4l. (Q(IsJ)sJ219,36) 
DO 43 I=1>31 

mene 495 (RK loJ), J=)]512) 
ОО 45 I=1s31 

READ 499 (RK(I9J) sJ=13924) 
ШІСІ I=1s31 

READ 49s(RK(IsJ)s J=25936) 
PRINT 8 


FORMAT (1H1) 

DO 1 151,31 

PRINT 31» (P(IsJ)sJ21518) 
PRINT 8 

DO 2 1=1%31 

PRINT 31» {P(1»J)»J=19»36) 
PRINT 8 

DO 3 1=]%31 

PRINT 139(RK(1»J)»J=1918 ) 
PRINT 8 

DO 4 1=1%31 

PRINT 13»{RK(I»J)»J=19s36) 
ШПЕМАТ (1Х%18Ғ5,2//) 

MS 200 

ERROR = 1.0 

DO 10 N= ІЗІҒ 

RAIL = OO 

ОО 29 1=2.30 

no 29 J22,35 

DO 15 K=1913 
P(10»K)=P(119K) 


P{10»14) = P(11s15) 
P(9914) = P(9915) 
P(8s14) : Р(8915) 
DO 19 K= 15919 
Р(ТҺЭК) - Р(8%К) 

P( 7920) = P(8,21) 
06520) - Р(6921) 
(5,20) = Р(6»21) 


P(5922)=P(6922) 
2(2523)-Р(6>23) 


Р(5»21) - Р(6921) 
P(5»24)=P(6524) 

el o— 275) = Р( 65326) 
Р(4»25) = P(4526) 
P( 325) = P( 326) 


Р(2 25) 70٦0 


99 





ai. u mr не! 


100 


ШУЛЫ) - Р(3526) 
Р( 2.27) = Р(3,27) 
Р(2,28) - Р(3,28) 
Р(2+29) = Р(3%30) 


P(1929)=P(1930) 
DO 18 K= 1517 

18 P( 309K) =P(295K) 
лао 18 ) = Р(29 4,17) 


0 18) = Р(29+17) 
768518) = P( 28517) 
PSIG), = P(26s17) 
OIG K=-19,26 
omer (279K) = P(26;K) 

ШІ 527) с Р(26%26) 
EE 21) = Р( 26.26) 
22527) = Р( 25.26) 
23.27) = Р( 23,26) 


ОО 17 К=27.36 
17 P(2&sK) = Р(23.К) 
CONST - 106000./P(15J) 
ieee ilsJ)) 51529951 
pul (RK(19J))52229352 
azo N=(-CONST *XQiIsJ)- (RK(I+1sJjJ)+RK(I1sJ))%*P(I+1sJ)+*#(FRK(1—]5 9 K 
1 (1, /)) Р(1-1».7) + (КК(Т»./+1) + КК (Т».))) ХР(То./+1) + (КК (Т».1-1) + ВК Та ру 2 


Ај 1) ) / (RK(I+1sJ)+RK(I-—1sJ)+RK(IsJj+1)+RK(IsJ-1)+&,* 
ОДА ТЈ) ) 
Ши ПЕККОК-АВ5Г СРЖЕРОТУЈ)) ) 26927527 
26 ҒАТ! - 1.0 
ШІР(ТІ».)) = PN 


29 CONTINUE 
BEINT 31» FAIL 
ШӘ (ҒАТ!) 33.33.10 
WW CONTINUE 
PRINT 8 
БОРО 30 І-1531 
30 PRINT 31» (Р(То.)) ».7< 118) 
٣۰٠٦٦٣۷۹۸۲٢ (1X,18F5.0//) 
ERINT 8 
DO 32 I=1931 
PREZERTNT 31> (P(1»2)3J=19»36) 
62 FORMAT (12Ғ6.0) 
49 FORMAT(12F6.2) 
41 FORMAT (18F4.0) 
44 FORMAT (18F4.3) 
END 
END POTMAP 





APPENDIX IV 


FUTURE PERFORMANCE PREDICTION 





102 
FUTURE PREDICTION PROGRAM 


PROGRAM PREDICT 
DIMENSION A, B, T, Q, DELP, P 2 
COMMON A, B, TEMP, PC; TC, T, Q, K, DELP 


READ Input for Z Calc. Constants 


READ DATA for Reservoir Properties 


READ Input for Q Funct. Tables 


READ DATA for Past Performance 
GAS, DIST., WATER | 


INC = CPI 






DO 100 TI TET 





READ 205 P(I) 









CONTINUE 


READ Input for future Performance 
GR, DR, WR 







Print out data for Reservoir 
Properties 








103 


Print out Past Production Data 
Calculate Z 


C = .00514590 (Temp) 


Bg, = (0) Ê 


“ 1858064) (AK) (Time) 
> POR)? 


СЕВ = (G) (Bg, ) 








104 


(WR) (TIME) 
(DR) (TIME) (CONV) 
(GR) (TIME) 


EKAY = GT + DT 





Calculate Water Influx 
| WE = (WIC) (WD) 


C)G - EKAY) 
(GRB + (WT) (BW) - WE) 


Calculate Reservoir Pressure 


РО = 


105 





Print Output 
Time, Pressure, Influx, Gas Produced 


STOP 


0. | Print: All wells 
have flooded out 
<0 
GE 
< O 


Print 1 Reservoir Depleted 
STOP 


106 
CALCULATION OF GAS DEVIATION FACTOR 


Subroutine Z CALC (S,Z) 


Dimension 





РҮ, = В, 
PY, = (8200. 
PY, = (B,)(3¥ - 1) 
РУ, = (В) (5%? - 31) 
Б 4 4 2 
РҮ. = (В.) (35% - 30% +0 
PY, = (B,) (637° - 70Y? * 15Y) 


(8) 
(82) (X) 
2 

(B4) X^ - 1) 
(B.)(5X^ - 3X) 

4 4 2 
(B.)(35X - 30x" + 3) 
(B4) (63% - 70X + 15%). 





RETURN 








107 


CALCULATION OF DIMENSIONLESS WATER INFLUX 
Subroutine QFUNCT (DELT, QT) 


TT = DELT 


— — — — — — —— DO 3 = кб, 1 


9 


QT 





ото о ТЕ ШЕШЕ ЕЕЕ 


B 


End 





108 


WATER INFLUX CONDITION 
Subroutine INFLUX (DELTD, WD) 


DIMENSION 









DELT = (TI)(DELTD) 
Call QFUNCT (DELT, QT) 









WI = (DELP ) (QT) 





ши — w и  “ “m Q Q ا س‎ a——  >—O— س‎ 


Es = 





CONTINUE 


RETURN 





109 
PRESSURE CALCULATION 
SUBROUTINE PRESS (Р, РОЙ, Р) 


DIMENSION | 
COMMON 















IF 
ABSF(E) - 1.0 


IF 
ABSF(E) - 0.1 


Zo 





со 


ESOS 


> 
=0 
<0 







ABSF (E) - . 


RETURN 


END 


110 








800 
101 


102 


200 


1000 


201 
1010 


202 
1020 


21013 


204 


205 
1030 


PROGRAM PREDICT 

DIMENSION A(696) 3B(6)5T(471)2Q(471) sDELP(999) »P(40) 
COMMON AsBsTEMPsPCoTCeTsQsKsDELP 

BRENT 800 

FORMAT(1H1) 

PRINT 101 

FORMAT(^8H RESERVOIR COMPUTER PROJECT---FUTURE PERFORMANCE>»//) 
PRINT 102 

FORMAT( 19H RESERVOIR DATA FOR>/) 

DATA INPUT FOR SUBROUTINE ZCALC CONSTANTS 

BON TODO J=1>6>1 

ШӘ ТОО 15156.1 

READ 200sA[19J) 

BORDUNSTCISX.F14,]]) 

CONTINUE 

ро 1010 151,6,1 

READ 201,B(1) 

FORMAT(F140.11) 

CONTINUE 

DATA INPUT FOR SUBROUTINE QFUNCT TABLE 

DO 1020 I=1>47193 

READ 2023T(I1)>0(I1)9T(I+1)90(1+1)>9T1(1+2)9Q(1+2) 
ЕОЮМАТ(Ғ12,2>Ғ12.23>5Ғ12.2>»Ғ12.,3»Ғ12.2»Ғ12.3) 
CONTINUE 

PANA INPUT FOR RESERVOIR PROPERTIES 


READ 
READ 
READ 
READ 
READ 
READ 
READ 
READ 
READ 
READ 
READ 
READ 
READ 
READ 
READ 
READ 
READ 
READ 
READ 
READ 
READ 


203sAIWEL 
203»СОҚУ 
203» ГТ 
203 „ОЕРТН 
203, ТЕМР 
203 1ہ‎ 
20390 
2039AH 
203956 
2039U 
203 AK 
9790 ۳ 
2005s. ٤۴ 
2039SWC 
203»РТ 
203 PC 
РТС 
203» ВИ 
203 „СРТ 
203» МІС 
20396 


FORMAT(F16.7) 


DATA 
READ 


INPUT FOR PAST PERFORMANCE 
204 »ОСА5 :ОТ5Т „МАТЕЕВ 


ҒОКМАТ(Е16.8»>4Х>Е16.8»4Х»Е16.8) 
INC=CPI 


DO 1030 I=1sINC,1 

READ 205,P(I) 

FORMAT(F10.23) 

CONTINUE 

DATA INPUT FOR FUTURE PERFORMANCE 











“amis a a ж. 
n + тош ыда» эт isl 


ин "LA 1, 


Eee un А Дон e Ú wa СЕ is 


"Wd М gy 
Oo mum wd эз 


A AMA % ЕУ Т 41 
“lr li ipa Fea Zn 
eaa а Је - шаа OED d 


(WP ate see کے‎ 


700 


300 
201 
50 
BID 3 
304 
305 
306 
pO 
308 
309 
310 
211 
12 
3155 


314 


316 
ү 


Sans 


103 
104 
Әже 
520 


221] 


READ 700sKK 

FORMAT(13) 

READ 2033 GR 

READ 203» DR 

READ 2033 WR 

DATA PRINT OUT FOR RESERVOIR PROPERTIES 

۲۴٢٢۲ و300‎ МС 

FORMAT(45H INCsCURRENT NUMBER OF PRODUCTION INCREMENTS= 914) 
PRINT 3015AIWEL 

FORMAT(48H MAX. NUMBER OF PRODUCING WELLS AT ANY ONE TIME=>F6+.1) 
PRINT 302,»BW 

FORMAT(50H BW»FORMATION VOLUME FACTOR FOR WATER IN BBLS/STB=>F6.3) 
PRINT 303, СОМУ 

FORMAT(44H CONVsSCF OF VAPOR PER BARREL OF DISTILLATE=>F7.1) 
PRINT 304sRI 

FORPAT(48H RIsRADIUS TO INITIAL EDGE OF RESERVOIR РЕЕ Т= о ٦| 
۰۱ 9) ٦۹ 

Про 1 Ззан ОЕРТН»ОЕРТН ОҒ RESERVOIR IN FEET=,F8.1) 

PRINT 306»ТЕМР 

FORMAT(50H TEMP,TEMPERATURE OF RESERVOIR IN DEGREES RANKINE-2,F7.41) 
ere 207 »PC 

FORMAT( 36H PC,PSEUDOCRITICAL PRESSURE IN А 11) 

PRINT 308, TC 

FORMAT( 50H TCsPSEUDOCRITICAL TEMPERATURE IN DEGREES RANKINE=5F6e1) 
PRINT 309,SGI 

FORMAT(40H SGIsINITIAL GAS SATURATION AS FRACTION=>F5.3) 

0 31090 

FORMAT( 24H O,POROSITY AS FRACTION=»F5.3) 

۱۱۱١۱ 4 ۲۶ 

FORMAT(43H AHsAVERAGE THICKNESS OF RESERVOIR ІМ ҒЕЕТ-»Ғ6.1) 
INT 312956 

FORMAT (44H SGs9SPECIFIC GRAVITY OF GAS RELATIVE TO AIR=,F7.4) 
IN 313,1) 

PORMAT(46H UsVISCOSITY OF RESERVOIR WATER IN CENTIPOISE=>5F663) 
PRINT 314,AK 

FORMAT(47H AKS,WATER PERMEABILITY OF RESERVOIR IN DARCIESz,F7.44) 
PRINT 315» TIME 

FORMAT(33H TIMEsTIME INTERVAL USED IN DAYS=»F7.2) 

INT 16 ٤٥ 

FORMAT(41H SGR,RESIDUAL GAS SATURATION AS FRACTION=,5F5.3) 

PRINT 317,SWC 

FORMAT(42H SWC$.;CONNATE WATER SATURATION AS FRACTIONz,F5.e3) 

PRESENT 318>WIc 

FORMAT(45H WIC»WATER INFLUX CONSTANT IN BBLS. PER PSIA=»F10.3»//) 
EW PRINT OUT FOR RESERVOIR PRODUCTION 

ERTNT 103 

PORMAT(55H RESERVOIR PRODUCTION FOR FUTURE PERFORMANCE PREDICTION) 
PRINT 104 

FORMAT(1H ) 

PRINT 319,GR 

FORMAT(41H AVERAGE DAILY DRY GAS PRODUCTION IN SCF=>F16+1) 

PRINT 320»DR 

FORMAT(44H AVERAGE DAILY DISTILLATE PRODUCTION IN STB=>F16.1) 
PRINT 321+WR 

FORMAT(39H AVERAGE DAILY WATER PRODUCTION IN STB=sFl6e1s//) 

MAIN BODY OF PROGRAM 





HEU 5 


322 


400 


40] 


106 


107 


402 


108 


یا ی ی۔ ری 
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S=P(1) 

БЕП ZCALC ( SoZ ) 

С-0.00514590%ТЕМР 

BGI=C#Z/5 
DELTD=1858064.0*AK*TIME/LU?FO*RI*#RI) 
GRB=G*BGI 

pe —=().0=-SGR-SWC)/SGI 

و٣۱۱‎ 600 ۴ 

PRINT 105 

BORMATLTISH INITIAL CONDITIONSs/) 

PRINT 322 sPI 

FORMAT( 26H INITIAL PRESSURE IN PSIA=»F8.2) 
223536 

eL OH INITIAL GAS IN PLACE IN SCF=>E1698) 
PRINT 400s»FLOOD 

FORMAT(47H MAXIMUM FLOODABLE RESERVOIR VOLUME IN BARRELSz,E16.8) 
PRINT 40152 

PORMAT(S86H INITIAL GAS COMPRESSIBILITY FACTOR=.E16.8377 )} 
ARENT 106 

١۷۸٣٢٦٢۳٦ »/»18H FUTURE PREDICTION»/)‏ ۰م 

PRINT 107 

FORMAT( 15H TIME PRESSURE :9X s2HWE216Xs3HGPWs8Xs/7HDELTA P) 
۳۳٣٣٣٣۲٦) 0٭‎ 0 

DO 2 و22‎ 1۱1 

ШЕСЕГТІ)-Р(І-І)-Р(Т) 

P1=P( INC) 

WT=WATER 

DT=DIST 

GT=GAS 

W=WR#TIME 

ШЕБЕСТІМЕЗЖСОМУ 

GDP=GRRTIME 

K=INC 

К=К+1 

LzK-1 

WT=W+WT 

DT=D+DT 

GT=GDP+GT 

EKAY=GT+DT 

DELP(K)20,0 

а ЕЕ TNFLUX(DELTD>WD) 

WE=WIC#WD 

POZ=C*(G-EKAY)/(GRB+WT*BW-WE ) 

CALL PRESS(P1»POZ»P2) 

DELP(K)=P1-P2 

PE ENT 402sLsP2eWEsEKAY sDELP (K) 

FORMAT( I15»2X»F8.2»2X»sE16.8»2X»E16.8»2X»F8.2) 
P1=P2 

IF(WE-FLOOD) 85454 

PRINT 108 

FORMAT(1H s/s27H ALL WELLS HAVE FLOODED OUT»/////) 
CONTINUE 

ЈРОКК—]) 6+6%»7 

CONTINUE 

GO TO 1 

IF(P1-14.7) 1031039 





10 
109 


(ОКО 3 
۱غ‎ 109 
FORMAT(1H 
CONTINUE 
Ger TO 5 
END 


As IS ERESERVOLRSBELBFTEDS////F) 
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SUBROUTINE ZCALC(Ss7Z) 

DIMENSION A(696) 2B(6) 9T(471) 9Q(471) sDELP(999) »P(40) 
DIMENSION PX(6)»PY(6) 

COMMON AsBsTEMPsPCsTCsT sQsK sDELP 
Y=((2.0*TEMP/TC)-4.0)/1.9 

ОО ЕВС) 

Ришар - ага | ку 

es) =6(3)%((3.04Y#Y 1-120) 

ШГ(0)-В(4)%( (5,0%Ү%Ү)-340) <Ү 

My (5)=B815)*( 1135.0" Y*Y )—-30.0) *Y* 73.0) 
map S)*C(l(623,05Y*Y)- TOQO) Sw 2 ET 5.9) *Y 
Х-((2.0%5/РС)-15.0)/14.8 

РХ(1)=В(1) 

PX(2)2B(2)*X 

Poems =o 3 jie (3.0*X*X)-1.0) 
PX(4)=B(4)*((5.0%*X*X)-30.0) *X 


565) 5(5)%((((35,;0*#X%#X)—30,0)%*%%xX)r5;0) 
ПРИ СОЈ КЕ Г 163 ХА РОС О УУ ТВО J xx 
APP=0.0 


DO 1 M=1>6931 

DO 1 N=1>6931 
APP=(A(N+»M)*PX(M)*PY(N))+APP 
Z=APP 

RETURN 

END 





SIBROUTTNE GFUNCTUBELTs OT) 

DIMENSION A(696) »B(6)9T(471) 3Q(471) sDELP (999) 5P(40)} 
COMMON AsBsTEMPsPCsTCsT sQsKsDELP 

TEMEEDELT 

ВО 3 1=1,471.1 

ТЕТО) 152.3 

И [—1)—Q(I-—2))Z(T(I=1)—T(1=—2)) 

ЕЕ Э ОТТЕ ТОА ССТ) ТИЕТ ИВО АСТ TCI) 
ИО [= 2)+(ттегт(122) ухројесттет стег у је ТТ тело а је 
6 ТО 4 

QT=Q(1) 

GO TO 4 

CONTINUE 

RETURN 

END 
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SUBROUTINE INFLUX(DELTD+»WD) 
DIMENSION A(696) 29B(6)2T(471) 2Q(471) sDELP(999) »P(40) 
COMMON AsBsTEMPsPCsTCoT sQsKsDELP 
М2= 0.0 

DO 1 J=1>K>1 

WJ إا‎ 

DELT=TI*DELTD 

Game QFUNCT (DELT sQ7T) 

N=K+1-J 

WI=DELP(N) *QT 

WD=WI +WD 

CONTINUE 

RETURN 

END 





ӘСЕР /ТІМЕ РКЕЗЗІРІ»РО2»Р2) 
ОІМЕМ5ІОМ А(6%»6)эЭВ(6)>Т(8471)>»0(471)>»ОЕГР(999)»Р(40) 
COMMON AsBsTEMPsPCsTCsT sQsKsDELP 
S25 1-5-0 
Grime ZCALC(SsZ) 
РО7Т=5/7 
9 POZ 
IF(ABSFIE)-5.0) 55552 
IF(E) 491753 
S=S-5.0 
GOTO 1 
S=S$+5.0 
COTTO 1 
ВЕБЕ Е) -1.0) 9.9.6 
ВЕНЕ) 9517.7 
5-5-1.0 
ЫШ ТО 1 
SzS41,0 
cio Or 1 
ШЕГАБ5Г(Е)-0.1) 13513510 
ШЕГЕ) 12.,17.511 
5-5-0.1 
БОР ТО 1 
Seo +) « | 
COSTO l 
ПЕ ГАБЭҒГ(Е)-0,02) 179175914 
ШИҮГЕ) 16%17515 
5=5-0.01 
СОТО 1 
5=5+0.01 
COTTO I 
Р2=5 
RETURN 
END 
END PREDICT 








€ DATA 


0.7071068 
1.224745 
0.7905695 
0.9354145 
0.265165 
0.293151 


С DATA 


INPUT FOR ИВО ти и ЕС CONSTANTS 


2. 1433504 
02083176184 
-0.021467042 
-0.00087140318 
0.0042846283 
00016505543 
02 ЕЕ 
-0.13403614 
0.066880961 
—O 0271785261 
0008517229] 
-0 00 2 1520929 
Ооо пас ри 
00050393604 
0.0050924798 
0201055023365 
->0 9 
Оса 69پ ّ پ""ٰ‎ 
-0.052184040 
0.044312146 
-0 5ھ‎ 
0 6ھ‎ 
О: О ге 
ص‎ )) ۶۶ OU 
On eO 
0-۔‎ ھ٦٢٣٢‎ 54 
0.019262143 
=0: 0 15 210600 
0.0042910089 
—0 . 00081 302526 
-0.0053095900 
О ۹0 
Ое 000580039 
-0.0060114017 
0.09275 то 


۹ 


٥۱۷٠۷٦٦٠۹٣٠٠٠٠۹٠ ٠٦ ١0٥ ۰ ۶۴ 


00000000 
00000010 
00000025 
00000050 
00000080 
00000200 


00000000 
00000404 
00000689 
00001020 
9 7ئ0 
00092447 


(0009506061 
7 5 
9 06 
00000060 
00000090 
00 О 


Ооо ۶۸۶ 
00000520 
О0о 0 758 
00001140 
00001469 
00003202 


00000005 
00000020 
00900040 
00000070 
00000100 
00000400 
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00000278 
00900606 
00000898 
00801251 
00001569 
00003893 





00000500 
00000800 
00001100 
00001400 
00001700 
00002000 
002 0 
00002600 
00002900 
00003200 
00003500 
00003800 
00004100 
00004400 
00004700 
00005000 
00005300 
00005600 
00005900 
00006200 
00006500 
00006800 
00007100 
00007400 
00007700 
00008000 
00008300 
00008600 
00008900 
00009200 
00009500 
00009800 
95010500 
00012000 
90013500 
00015000 
00016500 
00018000 
00019500 
00021000 
140022500 
00024000 
00025500 
00027000 
00028500 
00030000 
00031500 
00033000 
00034500 
00036000 
29037500 
00039000 
00040500 
00042000 
00043500 
00045000 


00004539 
nonn6314 
00007940 
00099461 
000109713 
000122379 
00013684 
0090150 13 
09 ڈ۶‎ 
0 0 
00018845 
00020080 
058 
00022500 
00023684 
00024855 
0090276020 
00027174 
00028314 
00029443 
00030565 
009951 79 
000327735 
0000 9 
90034974 
006036058 
Пишете 
0ئ0(‎ 7 
0003592 72 
00040331 
00041385 
00042433 
00044858 
00049968 
00054976 
09058285 
00064737 
09-512 
90074226 
70078886 
00083497 
000882062 
00092589 
000970601 
00101549 
98 
00 7 
00114738 
G 119083 
00123403 
00127699 
O01172 
0092792273 
00140453 
90144664 
00148856 


00000600 
00000900 
00001200 
Оо 
00001800 
ОООО СО 
00002400 
0 оО 
00003000 
00003300 
001903600 
00003900 
00004200 
00004500 
00004800 
0 ۵ 
00005400 
00005700 
00006000 
90006300 
00005600 
00006900 
90007200 
09007500 
00007800 
00098100 
99098400 
99008700 
22009030 
020 09399 
70٦ 
00009900 
00011000 
а 
00014000 
00 15505 
00017000 
00918500 
00 0992 
0021500 
00023000 
00024500 
00026000 
0092 7200 
000922900 
00030500 
000520006 
0003 35:00 
00035000 
90036590 
00038000 
00039500 
00941000 
00042500 
00044090 
00045320 


00905153 
70006869 
00008457 
00009949 
000112386 
00012778 
00014131 
00015450 
00016742 
00018011 
00012259 
00020488 
DOO E 
00022897 
00024076 
00025244 
00026406 
00027355 
00028691 
00029818 
09030937 
00032048 
cC 1 
00034247 
(00 5 
00036418 
00037494 
140038563 
00039626 
00040684 
00041735 
00042781 
00046574 
00051648 
۷٦ 
ОООО Я 
00066 36 
0 6 
یٰ0‎ 
00080428 
00085023 
018083575 
00094090 
00093 571 
ООО ОТО 
0107437 
ЕЛЕК 2 7 
971 106159 
00120526 
00124838 
9512 962.6 
80 | 
OTIS > 
00141859 
00146064 
00150249 


99000709 
00001000 
00001300 
00001600 
00001900 
00002200 
00002500 
00002800 
00003100 
00003400 
00003 707) 
00004000 
00004300 
90004600 
10004900 
00005200 
00005500 
00005800 
0000 07, 
00006400 
00006700 
00007000 
00007300 
00007600 
00007900 
00008200 
00008500 
00008800 
20009160 
00009400 
00009 780 
00010000 
00011500 
00013000 
00014500 
00016000 
0 0, 
ہ۶۶‎ 70 
30020500 
00 ,٤ 
00023500 
00025000 
00026500 
00028000 
00029 900 
000327800 
00032500 
00034000 
0 S500 
0003 7000 
000 20 
00040000 
00041500 
00043000 
00044500 
90046090 
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00005743 
00007411 
00008964 
00010^34 
000 ЕЕ 
09073232 
00014573 
00015882 
00017167 
00018429 
0 
00020894 
0002210 1 
81910. 2/9) 
909024466 
000295659 
(0 ۷۹۷ 
0002 Е 
00029068 
0003501 7 
0003153505 
00032417 
0009977 
00034611 
000585677 
00000. ٢ 
00٦ 
0-00 
)8 0 ۲ 
00041034 
00042084 
00043129 
00048277 
OQQO S oM 
00058265 
000551 1 
00067928 
000726 1 
0007/ ۷۶ 
00081965 
00086545 
00091084 
00095588 
O01000 7 
00104495 
00108904 
00113284 
O01176 E 
00121966 
001 26270 
00230550 
00134808 
00139045 
00143262 
00147461 
Шо Am 








00046500 
00048000 
00049500 
00052000 
00054000 
80 7000 
00059000 
00062900 
00064000 
00067000 
00069000 
00072000 
00074000 
00077000 
00079000 
00082000 
00084000 
00087000 
00089000 
00092000 
00094000 
00097000 
00099000 
02192000 
00104000 
00107000 
0109000 
7 ۱ 0۵( 
00114000 
№11 7000 
169119000 
Ши 1122000 
00124000 
29127000 
00129000 
89132000 
00134000 
00137000 
19159000 
00142000 
00144000 
00147000 
00149000 
qu 55000 
о 2500 
00170000 
0177500 
4135000 
20192500 
00200000 
10207500 
29215000 
mo 222500 
00230000 
2237500 
00245000 


00153029 
00157184 
019) ез 
3291658163 
10173639 
Оена 
0 6 
9195208 
00200542 
00208502 
00213784 
002217670 
00226904 
0023472] 
00239912 
00247668 
00222219 
DIE DO b 
00255629 
00273274 
002782353 
00285498 
OU 2 OSS 
00298543 
00303550 
00-11055 
2027180755 
1۷ہ‎ 
00328480 
00255906 
00340843 
00348230 
00353144 
00360496 
00255256 
00372704 
00 Z 2 
00384859 
00389705 
0056596959 
00401786 
00409013 
00413820 
00428196 
00446077 
00463863 
00481558 
00499167 
)0 ۹6 
00534145 
(0551522 
00568830 
0058686072 
(0603252 
00620272 
00637437 


00047900 
00048500 
00050000 
00052500 
0ٰ ٰ ٤٦ 
000 2 о 
00060000 
NONE 2S0 
00065000 
0006 7500 
00070000 
00072500 
00 0۵ 
00077500 
00080000 
00082500 
00908 2000 
2005 1208 
00090000 
09082500 
00095000 
یٰ0‎ 00 
90100000 
6 70 
0 0 
00107500 
DUO NOD 
0ٰ۵ 
0115000 
00117500 
0012005000 
00122300 
00125000 
0000 100 
OCT 20000 
00132500 
001723900 
DO 9 
00140000 
00142500 
00145000 
00147500 
009150080 
201927200 
00165000 
0> 
00180000 
00167500 
00192000 
00202500 
00210000 
DO ۹۳ 
00225000 
00232509 
00240000 
00247500 


00154416 
5 
05052698 
00169549 
ھ١‎ 05/2 577 
00183124 
001789352 
09196544 
00203201 
0020983275 
00216417 
00222980 
00229514 
00226020 
00242501 
00248957 
O2 ЗВ 
ОССЕ о 
7 ۹٦ 
00274545 
00280888 
0۹ ٦ 
00293514 
00299799 
06 
001" TA 
О 945 
00324760 
003202278 
00237142 
00343308 
00349460 
GOO 2002 3/54 
00301720 
ОШ ЕЛГӘ ДВ 
619217057214 
09390093 
00386070 
00292725 
0023223167 
00404197 
00410214 
00416220 
00434168 
00452016 
00469771 
00487437 
ооа осоо 
00922220 
00592925 
O ہہ‎ 50 
00574585 
09591805 
00608965 
00626066 
00643113 


00047500 
00049000 
0005710099 
00053009 
00056000 
00958000 
00061000 
00063000 
00066000 
00068000 
09071000 
00073000 
00076000 
00078000 
00081000 
00083000 
00086000 
00988000 
00091000 
O00930 0 
00096000 
00098000 
00 
00103000 
00106009 
00108000 
00111096 
00113000 
no116000 
009 . ٦ 
00 ٦ 
00123000 
ОО 2С 
00128000 
00131000 
0013350090 
00136000 
00138000 
00141000 
00143000 
00146000 
00148000 
ООЛ ӘНІН 
00160000 
00 
00175000 
00182500 
00190000 
00 60 
00205000 
0 ۶ ۳۲ 
00220000 
00227 300 
00233000 
00242500 
00250000 


2 


00 ۷ 
00159945 
90165444 
00170914 
00179069 
00184473 
00 23 
1م‎ 6 
00205854 
00211145 
00219046 
00224289 
0022418 
00237218 
00245086 
00250245 
00257 x 
002620753 
002707 > 
0027956815 
00283420 
00288473 
00296070 
00301053 
00308567 
0051 2 
00327082 
00326000 
00333883 
00338 о 
00345770 
00350688 
00358048 
00362942 
003702 W 
00375185 
00382432 
00387282 
00394543 
003 
90406606 
00411418 
00422214 
00440128 
00457945 
00475669 
00493307 
005708581 
0029 99 
0054 т 
00563068 
00580262 
UOT 
00614672 
0 ۲ 
00648781 











00255000 00660093 00260000 ۰ ۹ 
00270000 00 ۸۶ 00272009 0 DOT 
00285000 00727449 00290000 00738598 
00300000 +070 00395000 UEM 
20515000 00794042 00320000 0205075 
00230000 00827088 00225900 00838067 
00345000 00859974 00350000 (ОР Ра 
00360000 00892712 00365000 00903594 
89575000 06725509 00380000 00936144 
(990000 ۰م‎ ۷9 003999050 00968566 
00405000 00920708 00410000 و0‎ 5 
00420000 010223138 00425000 07033028 
00435000 01054409 00440000 O 1065052 
00450000 01086390 00455000 OOS O24 
0046 5000 0171718257 00470000 01128854 
00480000 01150012 00485000 01160574 
00495000 6 0 0 0 4 ۲ 
20520000 01234203 00530009 01255141 
20550000 а 3 00560000 0131709 
00580000 00309225 00590000 0۶ 
00610000 01421224 00620000 01441820 
00640000 01482912 00650000 01503408 
00670000 01544305 00680000 01564706 
00700000 01605418 00710000 01623729 
00730000 QUOD 00740000 01686490 
00760000 01726859 00 770000 01747002 
00790000 ШЕТӘ 7222 00800000 3156972785 
00820000 01847336 00830000 Ol ٤, 
00850000 01907243 00860000 0192716565 
00880000 01966943 00890000 OI EOS 
20910000 02026438 00920000 02046227 
00940000 02085744 00950090 02105473 
00970000 02144878 00980000 021764555 
01000000 02 203061 01250000 02688967 
01750000 02633908 02000000 04095800 
03000000 05899508 03500000 06780247 
05000000 09563099 06000000 11047299 
072300000 13531457 08000000 14350121 
10000000 17586284 1.590000 CO2 

€ МОША INPUT FOR RESERVOIR PROPERTIES 

2.0 

108,0 

7290. 

11450. 

690. 

„ВО 

.28 

64.4 

2975 

.30 

e050 

180. 


.30 


0026 5008 
00280000 
00295000 
00310000 
00325000 
00340000 
00352005 
00370000 
00385000 
00400009 
00415000 
00430000 
00445000 
00460000 
00475000 
00490000 
00510000 
00540000 
00570000 
00600000 
00630000 
00660000 
00690000 
00720000 
00750000 
00780000 
00810000 
00840000 
00870000 
00900000 
00930000 
00960000 
00990000 
01500000 
(0ْ ۵ 
04000000 
07000000 
09000000 
1500000 
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00682640 
00716280 
00749725 
00782992 
00816090 
00849028 
00881816 
00914459 
00946966 
00979344 
01011595 
01043724 
01075743 
01107646 
01139439 
01111125 
071242222 
01276024 
01338486 
01400593 
01462383 
01523872 
01585077 
01646011 
01706688 
017687120 
یی O01527‏ 
01887298 
01947065 
02006628 
02065996 
02125184 
02184216 
03164780 
05005720 
07650096 
۲ی 12 
۲ 11989 
25380090 





0.15 
8838. 
662.26 
387.48 
1.03 
7.0 
7.52757 
606898622000. 
eof 207540E+10 
8838. 
8764. 
5159 7 • 
8509. 
8364. 
8219, 
8100, 
9 
25000000. 
1250. 
30. 


1,54221800E+06 


8.67600000E+03 
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RESERVOIR DATA FOR 


INCsCURRENT NUMBER OF PRODUCTION INCREMENTS= 7 

MAX. NUMBER OF PRODUCING WELLS AT ANY ONE TIME= 2.0 

BWsFORMATION VOLUME FACTOR FOR WATER IN BBLS/STB= 1.030 

Ban ScE OF VAPOR PER BARREL OF DISTILLATE= 768.0 

RIsRADIUS TO INITIAL EDGE OF RESERVOIR IN FEET= 2190.0 

DEPTH»sDEPTH OF RESERVOIR IN FEET= 11450.0 

TEMP TEMPERATURE OF RESERVOIR IN DEGREES RANKINE= 690.0 

PCsPSEUDOCRITICAL PRESSURE IN PSIA= 66223 

TCsPSEUDOCRITICAL TEMPERATURE IN DEGREES RANKINE= 38725 

SGIsINITIAL GAS SATURATION AS FRACTION= .800 

OsPOROSITY AS FRACTION= „280 

AHs AVERAGE THICKNESS OF RESERVOIR IN FEET= 64.4 

SG>»SPECIFIC GRAVITY OF GAS RELATIVE TO AIR= „7975 

UsVISCOSITY OF RESERVOIR WATER IN CENTIPOISEs .300 

AKsWATER PERMEABILITY OF RESFRVOIR IN DARCTES= .0500 

TIMEsTIME INTERVAL USED IN DAYS= 180.00 

SGR»sRESIDUAL GAS SATURATION AS FRACTION= „300 

SWCsCONNATE WATER SATURATION AS FRACTION= .150 

WICsWATER INFLUX CONSTANT IN BBLS. PER PSIA= 1528 

RESERVOIR PRODUCTION FOR FUTURE PERFORMANCE PREDICTION 

AVERAGE DAILY DRY GAS PRODUCTION IN SCF= 25000000.0 

AVERAGE DAILY DISTILLATE PRODUCTION IN STB= 1250.0 

AVERAGE DAILY WATER PRODUCTION IN STB= 30-0 

INITIAL CONDITIONS 

INITIAL PRESSURE IN PSIA= 8838.00 

INITIAL GAS IN PLACE IN SCF= 6.06898622E 11 

MAXIMUM FLOODABLE RESERVOIR VOLUME IN BARRELS= 2.22914288E 08 

INITIAL GAS COMPRESSIBILITY FACTOR= 1.32981926E 00 

FUTURE PREDICTION 

TIME PRESSURE WE GPW DELTAR 
7 8009.66 3,01997662E 05 2.54650962Е 10 90.34 
8 7870.52 3.92764487E 05 000۔3‎ 6 9 139414 
EM | 7734. 90 4.98711813E 05 3.48106962E 10 135.62 
№ 7602648 6.17296089E 05 aca E о 122002 
01 7473.08 ГЕ ЗОЗОЕ 05 4.41562962Е 10 129,40 
Ro 7346.57 8.89432311E 05 4.88290962E 10 126.51 
13 722 1 0 6 6 5.35018962E 10 7(9 а 
№ 7101.81 1.20466079E 06 8 10 1228192 
15 6983.39 ПИО ЗВОЗПЕ 06 6.28474962Е 10 вело 
16 6867.53 1.55959584Е 06 6.75202962Е 10 115.86 
№ 6754.18 1.75105577Е 06 72193008 >: 10 1135 
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2585. 79 
3512.46 
3459.66 
3407.38 
32395601 
3304,23 
22 210543 2 
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VITA 


Ossian R. Butterfield was born at Leominster, Massa- 
chusetts, on May 12, 1927, the son of Ossian R. and Helen P. 
Butterfield. He attended public schools in Leominster, 
graduating from Leominster High School in 1944. He entered 
the U. S. Navy in June, 1944, was enrolled in the Navy V-12 
College Training Program and attended Williams College in 
Williamstown, Massachusetts. In November, 1945, he was 
transferred to the U.S. Naval Reserve Officers Training 
Course Unit at Brown University. In 1947, he was graduated 
from Brown with the degree of Bachelor of Science in Engi- 
neering. Upon graduation, he was commissioned an Ensign in 
the U.S. Navy Civil Engineering Corps. Duty assignments have 
included engineering duties at the Naval Proving Grounds, 
Dahlgren, Virginia; the Naval Air Station, Dallas, Texas; the 
Naval Stations in Argentia, Newfoundland, and in Bermuda, The 
Mest Indies; the Naval Air Station, Patuxent River, Maryland; 
the Naval Radio Station, Kamiseya, Japan; and the Naval Mine 
۲۰٠٠٢۹٣۲٢ School, Yorktown, Virginia. In 1955 and 1956 he 
attended Rensselaer Polytechnic Institute in Troy, New York, 
where he received the Bachelor of Civil Engineering degree in 
June 1956. He has served as Executive Officer and as Command- 
me Orficer of U.S. Naval Mobile Construction Battalion Six 


deployed to Rota, Spain, and other bases in the Atlantic Ocean 
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Area. His most recent duty assignment was in the Bureau of 
Yards and Docks, Navy Department, Washington, D.C., where he 
was Director of the Construction Division responsible for the 
management of projects engineered and constructed by the U.S. 
Navy for other governmental agencies. 

He was selected for postgraduate education in petroleum 
engineering and entered the Graduate School of The University 
of Texas in September, 1963. 

Butterfield was married to Ann Churchill of King George, 
٣۰٣٠٢٢ ہ,‎ on April 11, 1950. They have two children, Marcia. 
age 13, and Brian, age 12, 


Permanent address: 6600 Willamette Drive 
Austin, Texas 


This thesis was typed by Virginia F. Dailey. 
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